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Element  Pattern  of  an  Axial  Dipole  in  a  Cylindrical 

Phased  Array,  Part  1:  Theory 


1.  INTRODUCTION 

Circular  cylindrical  antenna  arrays  are  of  interest  for  radar  and  communica¬ 
tion  applications  because  of  the  uniformity  of  their  circumferential  radiation  char¬ 
acteristics.  Some  cylindrical  array  applications  require  a  high  degree  of  per¬ 
formance  antenna  characteristics  such  as  precis  beam  pointing,  accurate  angle 
tracking,  ultralow  sidelobes.  and  wide  bandwidth.  To  achieve  the  associated  tight 
array  illumination  tolerances,  mutual  coupling  between  the  radiating  elements 
must  be  accurately  accounted  for. 

Surveys  of  the  state-of-the-art  of  cylindrical  phased  arrays  and  their  design 

1  2 

considerations  are  given.  *  However,  most  of  the  cited  bibliographies  deal  with 
aperture  arrays.  On  the  other  hand,  a  dipole  element  is  an  attractive  choice  tor 
an  array  radiator  due  to  its  simplicity  of  manufacture,  reasonably  wide  bandwidth, 
and  polarization  purity. 

These  factors  motivated  the  two-phase  effort:  theoretical  and  experimental. 


(Received  for  publication  31  January  1984) 
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tennas.  A. A.  Oliner  and  G. H.  Knivtel,  Eds.,  Artech  House,  Dedham,  Mass. 

2.  Mailloux,  R.  J.  (1982)  Phased  array  theory  and  technology,  Proc.  IEEE. 

70 (No.  3). 


1 


In  the  drat  phase,  cylindrical,  stacked,  ring-antenna  arrays  of  uniformly  spaced 
axial  dipoles  were  investigated  to  establish  the  mutually  coupled  element  pattern 
characteristics.  A  detailed  account  of  this  investigation  is  provided  in  Part  1  of 
this  report  and  has  been  briefly  reported.  *  The  second  phase  of  the  study  con¬ 
sisted  of  experiments  augmenting  and  supporting  the  theoretical  phase  and  is  pre¬ 
sented  in  Part  2. 

In  the  theoretical  analysis,  a  modal  approach  *  a  employed  that  utilises  a  unit 
cell  method  similar  to  that  reported  for  an  array  of  aperture  elements  in  a  circu¬ 
lar  cylindrical  ground. *  The  modal  approach  uaa  chosen  for  several  reasons. 

&  ft 

Asymptotic  treatments  '  may  be  well  suited  to  modeling  arrays  of  arbitrary  con¬ 
formal  shapes  after  their  numerical  validity  has  been  ascertained  by  comparison 
with  exact  solutions  of  relevant  canonical  problems.  The  modal  technique  avoids 
the  uncertain  accuracy  of  the  asymptotic  approaches  and  has  the  following  advan¬ 
tages:  (1)  In  principle,  it  furnishes  an  exact  solution  subject  only  to  truncation  er¬ 
rors  that  can  be  estimated;  (2)  It  is  uniformly  valid  in  all  spatial  regions,  whereas 
asymptotic  formulations  employ  different  representations  in  various  angular  do¬ 
mains  end  require  transition  functions  to  connect  the  results  smoothly;  (3)  It  has 
additional  flexibility  in  that  it  provides  answers  for  both  rectangular  and  isosceles 
triangular  array  lattices,  the  latter  after  a  minor  modification;  and  (4)  It  yields 
both  the  realised  gain  and  element  pattern  phase. 


2.  ANALYSIS  FOR  RECTANGULAR  LATTICE 
2.1  Array  Model 

The  cylindrical  array  model  under  consideration,  shown  in  Figure  1,  consists 
of  an  infinite  number  of  equispaced  stacked  rings  of  axial  dipoles.  The  rings,  of 
radius  pQ.  contain  N  equispaced  identical  dipoles  each  located  coaxially  s  distance 
s  above  an  infinite,  perfectly  conducting,  circular  cylindrical  surface  of  radius  s. 


3.  Herper,  J.C. ,  Hessel,  A,,  Mandarino,  C.,  and  Tomasic,  B.  (1980)  Perform¬ 

ance  of  a  dipole  element  in  a  cylindrical  array  -  A  modal  approach,  IEEE 
Antennas  Propag.  Symp.  Dig. ,  University  of  Laval,  PQ.  Canada,  182-165. 

4.  Sure&u,  J.C. ,  and  Hessel,  A.  (1971)  Element  pattern  for  circular  arrays  of 

waveguide -fed  axial  slits  on  large  conducting  cylinders.  IEEE  Trans.  Anten- 
nas  Propag.,  AP-19,  64-76. 

5.  Eichmann,  G.,  and  Lee,  K.S,  (1978)  Analysis  of  conformal  scanning  dipole 

arrays.  Antennas  Propag.  Symp.  Dig, ,  University  of  Maryland.  Washing¬ 
ton,  D.C. 

6.  Indenbom,  M.  V. ,  and  Filippov,  V,  S.  (1978)  An  asymptotic  solution  of  the 

problem  of  mutual  coupling  between  the  radiators  of  s  convex  cylindrical  ar¬ 
ray  antenna.  Radio  Eng,  Electron  Phys. .  23(No.  1),  42-49. 
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Figure  1.  Geometry  of  the  Circular  Cylindrical  Array  of  Dipoles  in  a 
Rectangular  Lattice 


The  dipoles  have  a  length  L.  narrow  width  w,  and  negligible  radial  thickness. 

Their  central  gaps  of  length  h  are  excited  through  identical  matching  networks  via 
respective  transverse  electromagnetic  (TEM)  transmission  lines  fed  from  matched 
generators.  For  the  case  of  single -element  excitation,  all  but  one,  of  the  genera¬ 
tors  are  short-circuited.  A  typical  wedge  -shaped  unit  ceil,  of  axial  dimension  d 
and  width  b  Inp^lV  at  the  dipole  location,  is  seen  in  Figure  2,  which  shows  a 
dipole  without  its  feer  . 
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Figure  2.  Unit  Cell  Geometry  of  a  Cylindrical 
Dipole  Array 


2.2  Active  Dipole  Current 

The  performance  of  a  singly  excited  dipole  in  a  match  -  terminated  array  envi¬ 
ronment  is  determined  by  superposition  of  fields  generated  by  the  various  discrete, 
independent  azimuthal  phase  sequence  excitations,  as  well  as  by  those  of  the  con¬ 
tinuous  spectrum  type,  the  latter  corresponds s  to  uniform-amplitude,  progres¬ 
sive-phase  axial  beam  steering.  A  combination  of  one  of  each  kind  constitutes  an 
active  cylindrical  array  excitation. 

As  a  first  step  toward  analysis  of  a  free-excited  array,  one  considers  its 
forced  active  excitation.  The  generic  boundary  value  problem  in  this  case  reduces 
to  that  of  evaluation  of  fields  radiated  by  a  dipole  in  a  radial  unit  cell  With  phase* 
shift  wails  and  excited  by  a  prescribed  gap  voltage.  Such  analysis  leads  to  the  de¬ 
termination  of  the  unmatched,  active  dipole  gap-impedanoe  as  a  function  of  circum¬ 
ferential  and  axial  phasing,  frequency,  and  geometry.  With  this  information,  an 
appropriate  matching  network  may  be  configured  and  the  element  pattern  evalu¬ 
ated. 

To  begin  with,  one  observes  that,  for  w/^«  1,  the  usual  assumption  Of 
strictly  axial  dipole  currents  should  yield  a  good  approximation  to  the  desired  di¬ 
pole  performance.  Asa  consequence,  all  fields  may  he  derived  from  a  Single 
axial  component  of  a  vector  potential  Az.  The  treatment  is  facilitated  by  introduc¬ 
ing  a  unit  cell  Green's  function  G(r,  r';  v,  kzQ)  which  represents#  to  within  s  con¬ 
stant  factor  the  vector  potential  due  to  an  axial,  electric  point  current  element 
located  at  r  =  t  in  a  radial  unit  cell.  Here,  v (v  *  0,  1, . . . ,  N  -  1)  denotes  the 
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circumferential  Pth  phase  sequence  excitation  with  a  uniform  unity  amplitude  and 
a  p.iase  delay  once  around  the  array  of  2irv  rad.  The  argument  kz0  indicates  that 
the  progressive  phase  delay  between  the  neighboring  elements  in  any  axial  column 
of  the  active  array  is  k^d. 

By  linearity,  the  expression  for  the  vector  potential  in  a  unit  cell  due  to  the 
dipole  surface  current  density  K  may  be  written  in  the  form 


A,(c*'.kzo)  =  Mo/ G  (l.r>.kzo)  K(r')  dS*  (1) 

where  r  -  (p,  4>,  z)  and  S  is  the  dipole  surface.  Here,  G  satisfies 

V2  G  +  k2  G  =  -  S(r  -  r').  k  =  ~  (2a) 

and  is  subject  to  the  following  boundary  conditions: 

G(a,  0  .  z;  r-)  =  0  (2b) 

radiation  condition  for  p  •*  «  (2c) 

two  Floquet  conditions: 

G(p.^-.  z)  =  o_iw*  G (p.~|  ,zf  <2d> 

G(p.<£.-y)  =•  G •  (2e) 


For  e^w*  time  dependence,  the  various  field  components  are 

jto  /32A_  \ 

6'=~(-3^  +  k,A') 


k2p  3<t>  dz 


H 


1_  3A* 

Mo  dp 


(3b) 

(3c) 


From  Eqs.  (3a),  (3b),  and  (1),  it  is  evident,  in  view  of  Eq.  (2b),  that  Ez  and  E0 
vanish  at  p  -  a.  Also,  from  Eqs.  (1),  (2a),  and  (3c),,  one  finds  that 
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H*(p$.#.*)  -  -  KM)  . 


(4) 


Finally,  to  satisfy  Eqs.  (2d)  and  (2e).  G  is  expanded  in  ternris  of  a  complete,  ortho* 
normal  Floquet  basis 

$mn(0.z)  =  -p-  +  ’Sm*1  (5) 

v«a 

where  i/m  *  v  +  mN,  (  v  •  0,  1, . . . ,  N  -  1;  m  s  0,  ±1, . . . )  and  kEn  ■  kzQ  +  2rrn/d, 
(n  =  0,  H....);  Setting 


G(f,r>.klo)*  2  Bm„(p.p'Ji>.k10)  4>mn(^.*)  4>S,n  (<*>'.*') 


one  finds  from  Eq.  (2a)  that  the  radial  Green's  functions  gmn(P»  P')  obey 


1  d  da_„ 


P 


To  satisfy  Eq.  (2b).  one  requires  gmn(a»  p  ')  -  0,  and  to  insure  Eq.  (2c). 
gmn<P*  P")  must  have  the  form  of  an  outgoing  wave,  or  decay  as  p--.  With 
square  root  choice  Re  [kR]  *  Re  [(k2  -  k2^*  2]  >  0.  and  Im  [kR  J  ^0,  the  solu¬ 


tion  for  g_„  is 
mn 


Omn  «P*P')  = 


2| 


Z*'m(K"P<)  =  J*'m,K«p<)  “  H^(KhP<) 


and  ■  max  (min)  (P,  P').  The  and  H„J  denote  Bessel  functions  of  the 

first  kind  and  Hankel  functions  of  the  second  kind,  both  of  order  v  . 

m 

The  forced  dipole  excitation  implies  a  prescribed  gap  voltage  Vg  and  the  asso¬ 
ciated  gap  electric  field  E_  «  E „*rt,  such  that  V  «*Ih. 

g  g  U  g  g 

the  boundary  conditions  on  the  dipole  surface,  E  *  0  on  the  dipole  arms  and 
Ez  a  Eg  in  the  gap,  lead  via  Eqs.  (1)  and  (3a)  to  a  linear  integral  equation  of  the 
first  kind  for  the  unknown  dipole  current  density  (assumed  to  also  flow  through  the 
gap): 
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(8a) 


QK  =  f  6  (£.£')  K(£')  dS'  = 

s 


0  on  dipole  arms 
Eg  in  tha  gap 


where 


*  jo>Mo  32  . , 

Q  (£.£*)  =  “  *  k2)  Q(L£')  .  (8b) 

The  integral  Eq.  (8a)  is  solved  approximately  via  Galerkin's  procedure  by  first 
expanding  K(r')  in  the  form 

o 

K(t')=2  c  V  (0‘,z*)  (9) 

q=l  ^ 

in  terms  of  a  finite  set  of  linearly  independent  basis  functions  ip  ,  such  that 
#q(0-.  *0  *  0  for  \4'\>  w/2pQ.  |e*|>L/2.  <0^(0%  z')  =  tPq(z')  for  |  0  *|<  w/2pQ» 
\z'\  <  L/2  and  0^(0',  L/2)  *  0^(0',  -L/2)  =  0,  where  we  have  assumed  a  uniform 
current  density  across  the  dipole’s  narrow  dimension.  Using  Eq.  (r>a)  in  Eq.  (8a), 
one  has,  in  operator  form, 

o 

^  cq  §  =  PhPw  (10) 

where  p^  =  p^(z^ )  and  pw  =  Pw(0")  represent,  respectively,  a  unit  pulse  of  width 
h  centered  at  z'  =  0,  and  of  width  w  centered  at  4>'  -  0.  Relation  (10)  is  enforced 
by  taking  its  moments  with  each  of  the  |^.  This  procedure  yields  the  desired  set 
of  linear,  inhomogeneous  equations  for  the  unknown  set  of  c^,  that  is, 

a 

1  Cq  0?p.£¥q)  =  Eg  (?p.  PhPw)  (  p  -  1.  2 . P  =  Q).  (11) 

The  inner  (scalar)  product  is  (f,  g)  defined  by 
w/2  L/2 

(f.fl)  =  f  d<^  f  9(0.*)  d*  •  (12) 

-w/2po  -L/2 


Explicitly,  Eq.  (11)  reads 


flfoPo  S 

'  2  q  =  i 


^^(^nPo)  2t>m(knp0)  (N^p,  $mn 


=  l^p*Ph£w) 


(13) 
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with  f  ■  (V‘0>1/a  "  iao*<0*  Choosing 


tin  “  {*  + j),  q*1.2...,a  •*»*  1*1  * -~* 


(14*> 


one  find? 


W7Tt 


=  poV5aSvmC^ 


where 


%* 


t>4c) 


c«**' 


-  e^j 


tm) 


Furthermore, 


«"*Pb  h/* 


fp»  PhPw)  =  «W>jf  ~  Bp 


€J5*> 


where 


8_  »  si*  . 

p  2  wnfih 


«MW 


Substitution  of  Bqs.  (14)  and-  (15*  into  !<%.  (t5>  yields  the  foltoudng  system  «d  Wlm 
inhomogeneous  equations  for  the  determination  of  tile  unknown  dipole  Giurrent  ex¬ 
pansion  coefficients: 


*  Aw<V>  =  ®p 


(p  =  Ik  •  •  •  ►  F  *  Q)'- 


where 


8 


(16c) 


2.3  Active  Gap-Impedance 

Once  the  dipole  current  density  K(z)  has  been  arrived  at.  the  active  gap-im- 
pedance  is  determined  via 


(17) 


where  <Ig>  denotes  the  average  gap  current.  Substituting  Eq.  (10a)  into  Eq.  (17) 
using  Eq.  (16c).  one  has 


Za^-kto)  =  B90'.l4o>+iVi'.iclo)  = 


N7T»L*k 

Q 

4d  2  (SqW)  B. 


(18) 


2.4  Matching  Network  and  Reflection  Coefficient 

For  a  good  radiation  efficiency,  an  identical  network  is  incorporated  in  each 
TEM  feed-line  to  match  the  dipole  impedance  Zg  for  a  selected  pair  of  values  of  v 
and  kZQ.  This  lossless  reciprocal  two-port  network  is  characterized  by  its  scat¬ 
tering  matrix  £,  It  is  defined  between  the  reference  planes  on  the  input  and  output 
transmission  lines  with  characteristic  impedances  Zj  and  Z2  and  is  illustrated  in 
Figure  3.  Three  specific  forms  of  matching  networks  including  derivation  of  re¬ 
spective  S-parameters  are  discussed  in  Appendix  D. 

With  the  aid  of  Figure  3,  one  can  determine  the  active  reflection  coefficient  r 1 
at  the  input  to  the  matching  network  to  be  (see  Section  2,  Appendix  D) 


i-Sj, 

where  t  -  SnS22  -  S12S21i 


(19a) 
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(19b) 


and 


V9{if,k«>)x  2,(1/. + 

is  the  unmatched,  active  gap  reflection  coefficient.  In  Eq.  (19b),  fi2  mnd  *j  <,*“ 
note  the  propagation  constant  and  the  length  of  the  output  transmission  line. 

The  voltage  Vg  across  the  gap  is  given  by  Eq.  (Z>14)  as 


V*'  k*o> 


S>i[l  +  -i*gf tt  v 

1-SM  ‘ 


(90) 


W*: 


V 


•it 

•ia 

Si 

St 

h*— **<* 


Figure  3.  Equivalent  Network  for  a  Dipole 
Element  in  a  Cylindrical  Array 


2.3  Active  Array  Fields 

From  Eqs.  (3a).  (1).  (6a),  and  (7b).  one  finds  for  P  >  P0  the  axial  hdffipbtteffl 
of  the  electric  field  of  the  active  array 


o 

2*  eq$V$iL^Wia.  (91) 

e  - 1 

Utilizing  Eqs.  (14b),  (15d),  (19)  and  the  relations  K  a  2n/ct,  b-  c ip^,  Kq.  (91) 
can  be  cast  into  the  form 


6t  [v.  **,)  =  —  1  <»•  Hi*  | Wl  • 

IDO  Bl.ft1"  * 


(92i) 
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where  the  active  transmission  coefficient  ie 


T*  <►.*»>  •  ~  -*'£>  (£]f  «»„  1 


(22b) 


2.6  Element  Pattern 

In  view  of  the  identity 


•m  tir  v/d 

So.fiot  *  JL  2  J**  #  r  t-*«d»  dh  la  *0.1 . N-1) 

”  2irN  i*o  J  •  tt*o*i  t 

-rr/d  '  o.t 


(23a) 


where  dQs  and  denote  the  Kronecker  deltas,  the  expression  for  E'®'  due  to  the 


'Ot 


(e) 


singly  excited  (s  ■  0,  t  ■  0)  element  becomes 

2  £*•*♦ 


41T**.  X  , 
«  w/d 


Transforming  the  integration  variables  k#n  ■  k^Q  +  2irn/d  -k^,  one  finds 


m 

-I 


E^d)  = 


jn  s-i  m 

~  ^11  .'I'-* 

*"Po  1*0  m> -m 


i  ,  *W  _ 

7  +Tn 


-3s 


^  /  Tmo  (l^P) •'** *  dk* 


_E  +  **  „ 
d  +Tn 


(24) 


where  *  =  (k2  -  k2)  1^2  and  where  the  Floquet  relation  Tmn  (v,  kzQ)  ■  Tm0 
( v»  kzn)  has  been  invoked.  The  sum  over  the  index  n  in  Eq.  (24)  may  be  further 
converted  into  an  infinite  integral  so  that 


N-1  «* 


E«*'  “  *  4^57  v<«  £  J..**"*  f  T~  <,<«»  ■».  • 


(25) 


/a\ 

Using  the  first  term  of  the  large  •argument  asymptotic  expansion  for  (tf  p), 
the  far  field  becomes 
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«r«~  vM  .•*  x 


(Mi) 


lm(Kp.x)  ■/«  Tmo^^t)  dk,  • 


(Mb) 


The  integral  in  Kq.  (28b)  Is  evaluated  by  the  stationary  phase  method.  telling 
*  k  coa  u,  P  «  r  sin  0,  i  *  r  cos  4 ,  stars  0  is  mstmrsd  from  te  cylinder  axis* 
one  may  write 

Kp  +  R,*  *kroo«(u-0)  mi 


which  results  in  a  relevant  saddle  point  •  0 . 

•kconi.  The  first-order  stationary  phase  avail 
yields 


tit  the  integral  ^  (»if.  s) 


W*Mf  *  2  (vACost) 


From  Eqa.  (26a),  (26b),  and  (26),  one  has 

1ST  ~  ®  v<«  **■*» 


(26a) 


where 


*•*>*  Tsztsn  t|t 


(26b) 


it  terns  out  feat  the  lowest  order  term  of  Eg  is  O  ~  (1/r2),  so  feat  fee  Eg  con¬ 
tribution  to  the  far  field  may  be  omitted,  and,  consequently.  Eg  constftutes  fee 
entire  far  electric  field  due  to  a  singly  excited  dipole  element  in  s  mateh-termi - 
nated  array  environment,  The  realised  element  gain  pattern 


(Ml|* 


(Ste) 
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defined  with  respect  to  the  available  power  Pinc“lVincla  /Zj  may  now  be  written 

as 


Denoting  the  complex  field  element  pattern  by 

«§WW*)«-|4*r5  F(0.0) 

Eq.  (29a)  can  be  expressed  in  the  form 


(30b) 


(31a) 


(31b) 


3.  ANALYSIS  FOR  ISOSCELES  TRIANGULAR  LATTICE 

3.1  Array  Model 

The  dipoles  are  now  arranged  in  an  isosceles  triangular  lattice  aa  shown  in 

Figure  4.  This  lattice  may  be  generated  from  a  rectangular  one,  for  example,  by 

♦r 

displacing  every  other  ring  by  w/N  rad  in  the  angular  direction.  The  location  of 
the  array  elements  on  a  cylindrical  surface  of  radius  pQ  is  defined  by  the  skewed 
curvilinear  basic  vectors  (see  Figure  4) 

p^  *  sp|  +  tpj»  (s  *  0,  1, « > . ,  N  *  Is  t  1  0i  Al,  ,  , , ),  (32) 

where  t  indicates  the  ring  and  s  the  element  number  In  it.  The  center  of  the 
singly  excited  element  (a  ■  0,  t  *  0)  is  located  at  (Pq,  d  *  0,  s  •  0).  The  extent  of 
a  wedge-shaped  unit  cell  of  rectangular  cross  section  is  2ir/N  x  2d  as  indicated  in 
Figure  4. 

3.2  Active  Gap-Impedance  and  Reflection  Coefficient 

The  analysis  exploits  the  previously  obtained  results  for  a  rectangular  grid 
geometry.  In  particular,  the  expansion  of  the  unit  cell  Green's  function  as  given 
in  Eq.  (6a)  also  applies  here  with  the  proviso  that  the  double  sum  ir  iices  are  con¬ 
strained  to  (m  +  n)  of  even  parity  as  demonstrated  in  Appendix  A.  Consequently, 
the  relevant  expressions  for  active  impedances  and  reflection  coefficients  in  an 
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HrfiA) 


bm«p0 


Figure  4.  Geometry  of  the  Circular  Cylindrical  Array  of  dipoln  in  a 
Triangular  Lattlca 


1 


isosceles  triangular  lattice  can  be  simply  inferred  from  die  solution  for  the  cor¬ 
responding  rectangular  grid  array  by  restricting  the  summation  hi  Eq.  (8a)  to 
<m  >  n)  even.  With  this  proviso  implemented  in  Eq.  (18b),  the  system  of  linear 
equations  (18a)  is  also  valid  for  isosceles  triangular  lattices.  Thus*  with  this 
modification,  the  relations  of  Sections  2. 3  and  2. 4  apply  to  cylindrical  dipole  ar¬ 
rays  with  isosceles  triangular  lattices. 
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3*3  Element  Pattern 

A  tingle  (a  *  0.  t  ■  0)  array  element  excitation  in  an  ieoacelea  triangular  lat¬ 
tice  of  Figure  4  can  be  expressed  in  terme  of  an  excitation  vector  g  *  { #cQ  #t0} 

(a  *  0.  I.....  N  *  If  t*0,  ±lt ...).  In  Appendix  B,  Jse  following  identity  ia  dem- 
onatrated: 


4»  4»  * 


(33) 


Each  term  in  Eq.  (33)  under  the  aum  and  integral  eigna  repreaenta  a  unit  •ampli¬ 
tude  progreaaive-phaae  array  excitation*  2ir*>/N  rad  in  angular  direction  and 
ks0d  +  rad  in  the  skew  direction.  In  view  of  Eq.  (33).  the  electric  field  due 

to  the  singly  excited  dipole  element  in  a  triangular  lattice  environment  ia 


/  ett;  ►. <*, 


(34a) 


-e/d 


where,  from  Eq.  (22a) 
V*a 


E,(r? 


(34b) 


(m+n)< 


and  Tmn  is  given  by  Eq.  (22b).  Following  the  procedure  of  Section  2. 6.  one  can 
readily  show  that  the  expressions  for  element  pattern  and  far  field  are  identical  to 
those  of  Eq.  (31).  where  now 


F  (».*)  = 


J 


TTkp.sin  6 


■jf  2  £  T^.e.koorf)  1 

1  V  v«-e  mi>a 


(35) 


and  is  given  by  Eq.  (22b)  with  n  *  0  and  klQ  *  k  cos  $< 


4.  NUMERICAL  ANALYSIS 

Based  on  the  above  analysis,  a  Fortran  IV  computer  program  was  generated 
for  evaluation  of  the  phase  sequence  active  gap-impedances,  of  the  phase  sequence 
active  reflection  coefficients  at  the  input  to  the  matching  network,  and  of  the  ele¬ 
ment  pattern  for  both  rectangular  and  isosceles  triangular  lattice  configurations. 
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il  Rectangular  Lattice 

The  dipole  current  expansion  coefficients  c^  [see  Eq.  (9)}  were  computed 

from  Eq.  (16a).  where  Q  is  ihe  number  of  current  terms.  The  (Q  x  Q)  matrix  in 

Eq.  (16a),  with  the  elements  A  as  given  by  Eq.  (16b),  is  Hermltian,  that  is, 

A  *  A*  ,  and,  therefore,  only  Q(Q  +  l)/2  matrix  coefficients  need  be  evaluated, 
pq  qp 

Furthermore,  one  observes  that  the  A  and,  consequently,  c^.  as  well  as  the  gap 
impedances  Zg(v,  k#Q)  and  the  reflection  coefficients  /^(f,  kaQ),  are  even  func¬ 
tions  of  v.  This  result  follows,  for  instance, -for  A^,  by  first  performing  a 
change  of  variables  v  -  v  -(N  -  l)/2  for  N  odd  and  v-*  v  -  N/2  for  N  even.  As  a  re¬ 
sult.  r -  ^N.“(N-i)/2  for  N  odd  and  **  I  ?-?N/2  ,or  N  From 

Eq.  (16b),  it  is  now  obvious,  upon  replacing  y  by  -  a  and  simultaneously  m  by  -m, 
that  A  tv)  •  A  <-*). 

pq  pq 

The  Bessel  and  Hankel  functions  in  expressions  for  A^  were  evaluated  by  the 
usual  numerical  methods,  ^  each  applicable  in  its  own  range  of  validity  With  respect 
to  the  argument  and  order.  To  avoid  numerical  overflow  or  underflow  In  regions 
where  asymptotic  expansions  apply,  the  Bessel  and  Hankel  functions  were  not  cal¬ 
culated  separately  but,  instead,  an  asymptotic  expression  for  was  em¬ 

ployed  as  described  in  Appendix  E.  hi  this  fashion,  the  exponential  dependence  of 
the  individual  factors  has  been  suppressed. 

The  double  sum  in  Eq.  (16b)  was  evaluated  in  the  form  Xngft  (XmfmH).  With 
respect  to  the  rate  of  convergence  of  the  series  in  Eq.  (16b)»  the  following  obser¬ 
vations  are  relevant:  The  series  £Q  converges  as  1/n3.  Since  w/2pQ  «  1  ih  the 
factor  Sw  as  given  by  Eq.  (14c),  the  rate  of  convergence  of  £_  is  predominantly 
determined  by  the  large  m  behavior  of  m  which  is  l/jr|m|  Me  The  slew 

convergence  of  this  series  with  respect  to  the  angular  index  m  was  accelerated  as 
described  in  Appendix  C.  With  such  convergence  acceleration  and  with  the  nine- 
significant-digit  accuracy  of  the  Bessel  functions,  a  numerical  accuracy  to  five 
significant  figures  in  Zg(v,  k£g)  was  obtained  for  ka  •  120,  using  modal  indices  m 
and  n  between  ±10. 

The  set  of  equations  (16a)  was  solved  using  Gauss's  elimination  method.  It 
should  also  be  mentioned  that  for  0  •  n/2  (ki  ?  =  0),  that  is,  in  the  principal  H- 
plane  cut  of  the  element  pattern,  the  axial  dipole  current  density  i®  reflection  - 
symmetric  with  respect  to  the  dipole  center.  For  an  arbitrary  phasing  k^g,  that 
is,  for  a  conical  cut  with  6  f  n/2,  odd  terms  also  appear  in  the  series  Eq,  (10a). 
However,  the  resulting  asymmetry  is  relatively  minor,  since  the  magnitudes  of 


7. 


Abramowitz,  M.,  and  Stegun,  I,  (1964)  Handbook  of  Mathematical  Functions. 
\J ;  S.  Dept.  Commerce,  N*t.  Bur.  Standi 
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even  (q)  dipole  current  terms  are  much  smaller  than  the  odd  (q)  ones.  This  fea¬ 
ture  is  due  to  the  constant  gap  E  -field  assumption. 

Using  the  known  dipole  current  coefficients,  the  active  gap-impedances  Z 
( v,  Hzq)  *  Rgfi',  +  J^g^»  kzo^  were  calculated  from  Eq.  (18).  Figure  5 
shows  the  "broadside  scan"  gap-resistance  Rg(0,  0)  versus  the  number  of  dipole 
current  terms  for  three  dipole  lengths.  L  =  0.  4A,  L  *  0.  5 A,  and  L  s  0.  6A.  The 
respective  reactances  Xg(0,  0)  are  seen  in  Figure  6.  Since,  in  this  case,  the  even 
dipole  current  terms  are  not  excited,  they  are  not  shown  in  Figures  5  and  6.  One 
observes  that  for  short  dipole  lengths,  a  single  current  term  yields  a  good  approx¬ 
imation  to  the  dipole  impedance.  It  was  found  that  the  curves  Z^(v>  0,  >  0) 

versus  Q  exhibit  a  similar  behavior  as  Z  (  v  =  0,  k  =  0)  shown  above. 

go 

The  relative  convergence  phenomenon0  was  observed  by  monitoring  the  behav¬ 
ior  of  current  expansion  coefficients  and  active  gap-impedance  values  as  a  func¬ 
tion  of  the  number  of  unit  cell  modes  n _ and  the  number  of  current  terms  Q.  It 

max 

was  found  that,  for  a  stable  solution,  it  is  necessary  to  impose  the  condition 
Q  <  2Lnmax/d*  For  example,  if  «max  *  15  and  d  =  0. 7,  one  requires  Q  <  17  if 
L  *  0.4  A,  Q  <  21  if  L  *  0.  5A.  and  Q  <  25  if  L  =  0. 6A.  If  Q  increases  beyond 
these  values,  the  numerical  values  for  tend  to  become  progressively  less  ac¬ 
curate.  as  can  be  seen  from  Figures  5  and  6.  This  condition  on  Q  insures  an  ade¬ 
quate  resolution  of  the  dipole  current  distribution  in  terms  of  the  unit  ceil  Floquet 
mode  basis.  Therefore,  the  highest  axial  spatial  frequency  of  the  truncated  unit 
cell  mode  basis  27rnmax/d  should  be  greater  than  the  highest  spatial  frequency 
Qtt/L  of  the  dipole  current  basis.  If  a  high  degree  of  accuracy  in  Z^  is  desired, 

that  is.  when  more  current  terms  are  needed,  one  must  also  increase  n _ ...  and 

max 


mmax*  respectively.  The  reason  for  the  increase  of  mmax  along  with  n 


max 


is 


the  behavior  of  the  Bessel's  functions  in  Eq.  (16b).  Namely,  with  increasing 

(2) 

n„,  .  *he  arguments  of  the  Z«  i*  :  a,  pn)  and  of  HV, '  (#„  Pn)  increase.  There- 

max  m  n  u  ^mnu 

fore,  to  achieve  a  sufficient  accuracy  in  A^,  it  is  necessary  to  increase  past 
the  transition  region,  that  is,  vm  >  |  *  n  pQ  | ,  where  the  convergence  of  the  series 
may  be  accelerated  as  described  in  Appendix  C. 

The  active  reflection  coefficients  r^(v,  kzQ)  at  the  input  port  of  the  matching 
network  were  computed  from  Eq.  (19a).  Two  simple  matching  networks  have  been 
considered:  (a)  series  and  (b)  parallel,  both  shown  in  Figure  7.  In  each  case, 

4 2  =  0  and  network  parameters  have  been  chosen  so  as  to  match  the  active  cylin¬ 
drical  array  at  v  -  0  and  kzQ  =  0  (see  Figure  3).  In  Figure  7,  the  transformer 


8.  Lee,  S.  W. ,  Jones,  W.  K.,  and  Campbell,  J.  J.  (1970)  Convergence  of  numeri¬ 
cal  solution  of  iris-type  discontinuity  problems.  Antennas  Propag.  Symp. 
Dig. ,  Ohio  State  University,  Columbus. 
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GAP  REACTANCE.  X.  (0,0)  (O)  SAP  RESISTANCE,  R„(0,0)  Ol) 


Figure  5.  Gap-Resistance  Rg<0,  0)  Versus 
Number  of  Dipole  Current  Terms  (b fk  •  0. 6, 
d/X  »  0.7.  s/X  «  0. 25,  w/x  ■  04  05.  h/X 
■  0.01) 


Figures.  Gap-Reactance  Xg(0,  0)  Versus 
Number  of  Dipole  Current  Terms  (b/X  *  0.  6, 
d/X  =  0.7,  s/X  =  0.25,  w/X  =  0.05,  h/X 
-  0.01) 
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Figure  7.  H-Plane  Voltage  Element  Gain 
Pattern  for  Cylindrical  and  Reference 
Planar  Arrays  and  Matching  Networks. 

(a)  Series,  (b)  Parallel  (b/X  *  0.  6,  d/X 
=  0. 7.  ka  =  120) 

ratio  nQ  ■  (Rg(0,  OJ/Zj)1^2  in  case  (a),  and  nQ  a  (Yj/Gg(0 ,  O)1^2  in  case  (b), 
where  Yg(v .  kz0>  “  llZg<V,  k2())  *  Gg(i/,  kzQ)  +  jBg(i/,  kzQ). 

The  element  patterns  of  the  excited  element  located  at  (pQ,  6  ■  90° ,  ^  ■  0°  ) 
were  calculated  in  the  spherical  coordinate  system  (r,  6,  4>  ).  having  its  z-axis 
coincident  with  that  of  the  cylinder  according  to  the  relations  in  Eqs.  (31)  and 
(29b). 


4.2  Triangular  Lattice 

The  Steps  in  Section  4. 1  apply  also  to  a  triangular  lattice.  The  dipole  current 
density  expansion  coefficients  sure  again  computed  from  Eqs.  (16a)  and  (16b)  but 
with  (m  +  n)  even.  The  convergence  of  the  series  in  Eq.  (16b)  was  accelerated  as 
described  in  Appendix  C.  Because  the  axial  unit  cell  dimension  is  2d  (see  Figure 
4),  the  relative  convergence  condition  is  Q  <  Ln  /d.  The  element  patterns 
were  calculated  from  Eq.  (31)  using  Eq.  (35). 

The  program  was  run  on  a  CDC  6600  computer  at  the  Courant  Institute  of 
Mathematical  Science,  New  York  University,  N.Y.,  and  at  the  AFGL  Computer 
Center,  Hanscom  AFB,  Mass.  All  calculations  were  carried  out  in  siugle  preci¬ 
sion  (15  significant  digits  on  a  CDC  6600  computer)  complex  arithmetic.  The  ex¬ 
ecution  time  for  a  180-point  element  pattern  is  less  than  one  minute.  The  memory 
required  for  this  program  is  150  kbytes. 


5.  NUMERICAL  RESULTS  AND  DISCUSSION 


The  amplitude  and  phase  of  the  element  patterns  were  computed  for  represent 
ative  values  of  array  parameters,  and  the  results  are  grouped  to  exhibit  the  sig¬ 
nificant  trends.  The  numerical  results  stress  the  element  performance  aspects 
relevant  to  cylindrical  array  design. 

In  order  to  maximize  the  broadside  element  gain,  a  matching  network  appro¬ 
priate  to  in-phase  excitation  of  all  dipoles  was  employed  throughout.  Except  when 
explicitly  stated,  the  following  values  of  dipole  length,  width,  gap  size,  and  dis¬ 
tance  to  the  cylindrical  ground  were  employed,  respectively,  L/A  *  0.5,  w/A. 

=  0.05,  h/A  *  0.01,  and  s/A  *  0.25.  In  all  cases,  field  amplitude  (voltage)  ele- 

2  1/2 

ment  patterns  were  normalized  to  the  unit  cell  gain  (4wbd/A  )  ,  and  10  dipole 

current  terms  were  used  in  Galerkin's  procedure.  Most  of  the  numerical  results 
were  generated  for  rectangular  grid  arrays  in  conjunction  with  the  design  and  per¬ 
formance  evaluation  of  the  cylindrical  array  described  in  Part  2  of  this  report. 
They  are  presented  in  Figures  7  to  20;  those  for  the  triangular  lattice  appear  in 
Figures  21  and  22.  Each  is  followed  by  a  detailed  discussion. 


Figure  8.  H- Plane  Voltage  Element  Gain 
Pattern  (d/A  3  0.  7,  ka  =  120).  Parameter: 
azimuth  spacing  b/A  *  0.5,  0.  6,  0.  7 
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5.1  Rectangular  Grid  Array 

5. 1. 1  INFLUENCE  OF  ARRAY  GRID  PARAMETERS 

5. 1. 1. 1  Element  Amplitude  Pattern  -  Forward  Region 

Figure  7  shows  a  comparison  of  the  principal  H -plane  cut  (  $  »  90“ )  of  the 
voltage  element  pattern  for  ka  =  120,  b /A  *  0.  6,  d/A  *  0.  7  with  that  of  the  corre¬ 
sponding  planar  array.  It  is  observed  that  both  patterns  are  similar  and  exhibit  a 
substantial  drop-off  near  ^  ■  42“  .  In  the  planar  array  case,  this  drop-off  is 
caused  by  an  end  fire  grating  lobe  condition  (EGL)  for  b /A  *0.6  and,  in  the  cylin¬ 
drical  array,  by  its  quasiplanar  counterpart.  As  expected,  the  planar  array  drop¬ 
off  is  steeper  (since,  in  view  of  the  curvature  in  the  cylindrical  array,  fewer  ele¬ 
ments  participate  in  the  EGL  effect). 

Figure  8  illustrates  the  b/X  dependence  of  the  9  •  90“ cut  for  ka  *  120  and  an 
axial  spacing  of  0.  7 A.  It  is  found  that,  as  in  planar  arrays,  an  increase  of  b/X 
causes  a  predictable  progressive  narrowing  of  the  element  pattern. 

The  curves  of  Figure  8  also  exhibit  a  ripple  in  the  broadside  region,  whose 
amplitude  diminishes  with  tightening  of  azimuthal  spacing  and  becomes  negligible 
for  b/A  =  O.S.  On  the  other  hand,  for  b/X  =  0.  7,  one  finds  a  peak-to-peak  ripple 
amplitude  of  about  0. 5  dB.  In  a  low  sidelobe  cylindrical  (or  curved)  array  design, 
it  may  be  necessary  to  reduce  the  ripple  amplitude  and  consequently  employ  a 
tight  azimuthal  spacing.  More  on  this  subject  will  be  said  in  Fart  2  of  this  report. 
The  ripple  is  not  an  edge  effect,  as  no  aperiodic  edge  discontinuities  exist  in  an¬ 
gularly  and  axially  periodic  cylindrical  arrays,  but  is  due  to  the  interference  of 
the  direct  single  element  radiation  (with  planar  element  pattern)  with  the  grating 

lobes  of  respective  quasiplanar  subarrays  excited  by  the  guided  creeping  wave 

1  4 

phase  gradient,  as  discussed  *  for  the  case  of  cylindrical  arrays  of  aperture  ele¬ 
ments.  Since  the  ripple  is  a  result  of  a  combination  of  the  angular  periodicity  of 

1  4 

the  array  and  its  curvature,  *  it  is  also  expected  to  appear,  and  is,  indeed, 
found  in  curved  periodic  dipole  arrays. 

To  exhibit  the  element  pattern  features  as  a  function  of  the  elevation  angle. 
Figure  9  shows  a  number  of  constant  9  (conical)  cuts  for  ka  =  120,  b/A  *  0.  6,  and 
d/A  *  0.  7.  Several  trends  are  observed  with  decreasing  values  of  6  measured 
from  the  cylinder  axis:  The  EGL  drop-off  broadens  and  shifts  toward  larger  val¬ 
ues  of  the  E-plane  gain  level  is  progressively  reduced,  and  the  ripple  ampli¬ 
tude  decreases.  The  first  two  of  these  features  may  be  explained  in  quasiplanar 
terms  and  are  predictable  from  the  local  planar  grating  lobe  diagram  adapted  for 
the  cylindrical  array  geometry  as  shown  in  Figure  9.  In  this  figure,  uQ0  -  v! ka 
=  sin  9q  sin  and  Uqq  -  k^/ k  *  cos  0q.  As  usual,  the  (m  *  0,  n  =  0)  reciprocal 

lattice  point  specifies  the  main  beam  direction;  all  other  points  (m,  n)  correspond 

2  2 

to  the  grating  lobe  location  in  (u,  v  )  space.  The  (m,  n)  lobes  with  u  +  v  <  1 
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are  visible.  Th«*  amount  cf  shift  of  SQL  dr«n-ofr  toward  larger  values  of  with 
decreasing  9  may  be  predicted  from  the  geometry  of  Figure  9,  which  leads  to  the 
relation:  * 


sin  4>  -  — t  1-  (cos  9  +  n  -  m-— -  1 

sql  sind  Iffll  1  d  /  b 


(m  *  0*  •  •  • )»  (n  *  0,  ill  •  •  •  )• 


(36) 


It  is  seen  from  Eq.  (36)  and  Figure  9  that  in  the  H-plane  (  9  •  9CT  ),  the  grating 
lobe  (m  *  -1,  n  «  0)  enters  into  real  spaoe  at  »  42*.  As  9  decreases*  the 

corresponding  shifts  toward  larger  values  of  causing,  as  already  men¬ 

tioned*  a  widening  of  the  element  pattern*  The  potential  EGL  drop-off  correspond¬ 
ing  to  the  axial  grating  lobe  (m  =  0,  n  *  -1)  for  d/A  *  0. 7*  as  predicted  from 
Eq.  (38)  is  not  noticeable  at  9  *  60,  *  25. 4*  The  tapering  off  of  the  EGL 

effect  with  decreasing  0  and  the  associated  diminishing  of  the  ripple  amplitude  in 
Figure  9  are  attributed  to  the  reduction  of  mutual  coupling  in  axial  dipole  arrays 
for  directions  off  the  principal  H-plane,  and  particularly  In  the  near  E-plane  re¬ 
gion. 

Figure  10  shows  three  E-plane  cuts  for  axial  spacings  d/A  *  0. 5*  0. 6,  0.  7, 
with  b/A  »  0. 6  and  ka  *  120.  It  was  verified  (down  to  ka  «  30)  that  the  axial  gain 
roll-off  with  6  matches  that  ih  appropriate  planar  arrays.  Although  an  axial  spac¬ 
ing  of  0.  7  A  admits  only  a  limited  (to  about  23°  off  broadside)  elevation  scan  of  an 
active  array,  element  patterns  for  large  values  of  9  are  useful  for  sidelobe  level 
evaluation. 

Figure  il  exhibits  the  dependence  on  cylinder  radius  (ka)  of  the  EGL  drop-off 
and  of  the  ripple  in  the  9  «  90*  cut.  It  is  seen  that,  with  decreasing  values  of  ka, 
the  EGL  drop-off  slope  diminishes.  Simultaneously,  the  ripple  frequency  is  re¬ 
duced  and  its  amplitude  increases.  The  tapering  off  of  the  EGL  effect  is  to  be  ex¬ 
pected,  since,  with  decreasing  ka,  fewer  elements  participate  in  its  formation. 

The  reduction  (with  decreasing  ka)  of  the  ripple  frequency  is  due  to  reduction  of 
the  creeping- wave  grating  lobe  phase  delay  variation  per  unit  increment  in  the  ob¬ 
servation  angle.  Finally,  the  decrease  of  ripple  amplitude  with  increasing  ka  is 
expected  due  to  the  increased  creeping  wave  loss.  In  the  planar  limit,  no  ripple 
is  observed. 


*This  relation  can  also  be  obtained  from  the  condition  of  transition  region  bound¬ 
ary  n^m  (#fnV»  that  ls*  I  ‘'m1  =  *nP0’  wlth  V*0P0  =  8in  *0  8ln  *0'  and 
kzn^k  =  cds  *  nX/d. 


22 


AZIMUTH  ANOLM  tt»0) 


Figure  9.  Voltage  Element  Gain  Pattern  and 
Corresponding  Local  Grating  Lobe  Diagram 
(b/X  "  0. 0.  d7X  *  0. 7,  ka  *  120).  Parameter: 
elevation  angle  0*  3(f,  40*,....  90* 
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Figure  10.  E-Plane  Voltage 
Element  Gain  Pattern  (b/X 
=  0. 6,  ka  s  120).  Parameter: 
axial  spacing  d/X  =  0. 5,  0.  6, 
0.7 


5. 1. 1.2  Element  Amplitude  Pattern-Shadow  Region 

Knowledge  of  the  shadow  region  slope  in  important  in  the  evaluation  of  far 

sidelobes.  Figure  12  illustrates  the  dependence  on  cylinder  radius  of  the  shadow 

region  slope,  in  the  0  «  9(P  cut.  It  is  seen  that  the  slope  is  linear  on  a  dB  scale 

and  is  found  to  increase  as  (ka)*^3.  The  linearity  of  the  slope  indicates  that  the 

angular  pattern  in  the  shadow  region  is  primarily  due  to  a  single  creeping  wave 

1/3 

contribution  with  an  angular  attenuation  constant  proportional  to  (ka)  .  The  rip¬ 
ple  in  the  ^  »  180*  region,  similar  to  that  found  in  the  case  of  au  unperforated  cyl¬ 
inder,  is  a  result  of  the  interference  of  two  creeping  waves  traveling  in  opposite 
directions  around  the  cylinder.  The  resolution  of  this  ripple  at  the  very  low  pat¬ 
tern  levels  attests  to  the  high  numerical  accuracy  of  the  theory  and  numerical 
techniques. 

Figure  13  presents  the  0  dependence  of  the  shadow  region  slope.  It  is  found 
that  in  the  range  of  0  shown,  the  slope  varies  as  sin1^3  0  .  The  overall  depend  ¬ 
ence  of  the  shadow  region  slope  on  (ka  sin  0)*^3  is  consistent  with  other  asymp- 

1  4 

totic  developments.  In  all  cases,  the  computed  shadow  region  slope  was  in¬ 
sensitive  to  axial  spacing. 

Figure  14  exhibits  the  dependence  of  this  slope  on  azimuthal  spacing,  for 
d/X  =  0. 7,  ka  *  120,  6  •  9(f  .  One  observes  (see  also  Figure  8)  that  the  patterns 
intersect;  while  the  element  pattern  for  b/X  =  0. 5  is  wider  in  the  broadside  region, 
its  shadow  region  slope  exceeds  that  for  b/X  »  0,7  so  that  its  shadow  region  ele¬ 
ment  gain  falls  below  that  for  b/X  «  0. 7.  This  result  is  not  unexpected,  because, 
fr  smaller  values  of  b/X,  the  creeping  wave  generating  the  shadow  region  H -plane 
pattern  encounters  a  larger  density  of  match-terminated  elements.  The  minor 
■jpendence  of  the  E -plane  pattern  on  b/X  is  not  shown. 

5. 1. 1.3  Element  Pattern  Phase 

Design  of  a  cylindrical  anr y  requires  knowledge  of  the  element  phase  center 
location.  In  Appendix  F,  we  describe  element  phase  pattern  reference  point  trans¬ 
form  ation  from  array  axis  to  the  point  ( Pph>  0  =  90° ,  0  =  or).  Namely,  it  was  found 
that  the  flattest  element  pattern  phase  in  the  broadside  region  is  obtained  for  the 
ph__,e  reference  location  approximately  at  ( p^  =  a  +  0. 8  s,  9  ■  90* ,  ^  *  0“ ),  that 
is,  0. 2  s  below  the  dipole  center.  Figure  15  exhibits  the  element  pattern  phase  in 
the  9  *  90*  cut  for  the  phase  reference  location  at  p^,  its  dependence  on  b/X,  for 
d/X  =  0. 7.  One  observes  that  the  phase  varies  only  by  a  few  degrees  up  to  the 
EGL  drop-off.  Thus,  the  principal  H-plane  EGL  position  essentially  determines 
the  limit  of  usefulness  of  the  element,  both  in  amplitude  and  phase.  Figure  16  pre¬ 
sents  the  E -plane  element  phase  pattern  versus  axial  spacing  for  the  choice  of 
parameters  of  Figure  10. 
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W*u»  19.  Element  Qain  Pattern  MB)  (b fk  ■  0.  f, 
d/A  *0.7,  ka  »  110).  parameter!  elevation  angle 
•  ■  1<P »  5<P  *  8 (f 


AZIMUTH  ANQLK.  AtOCQ) 


Figure  14.  H- Plane  Element  Gain  Pattern  (dB) 
(d/X  -  0, 7,  ka  «  120).  Parameter!  acimuth 
spacing  b/X  ■  0. 5,  0. 6,  0. 7 
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Figure  IS.  H- Plane  Element  Phase  Pattern  (d/1 
•  0. 7,  ka  ■  ISO).  Parameter:  aximuth  spacing 
b/X  «  0.8,  0.6,  0.7 


Figure  16.  E-Plane  Element  Phase  Pattern  (b/X 
■  0. 6,  ka  *  120).  Parameter:  axial  spacing  d/X 
*  0.5,  0.6,  0.7 
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5. 1.2  INFLUENCE  OF  DIPOLE  DIMENSIONS 

The  influence  of  the  dipole  length  L  and  of  the  distance  to  ground  a  on  the  cyl¬ 
indrical  array  element  pattern  was  investigated  in  the  9  *  9CP  cut  for  b  /A  *  0, 6 
and  d /A  •  0.  7.  For  each  selection  of  an  (L,  a)  combination,  the  array  was 
matched  for  in-phaae  excitation.  The  variation  of  the  element  gain  pattern  with  L 
and  (or)  a,  for  L/A  •  0.4,  0.8,  0.6,  and  0.2  C  a  Ik  4-  0.3  waa  found  inaignificant. 
On  the  other  hand,  the  influence  of  these  parameter*  on  the  gap  impedance  waa 
obaerved  (aee  Figure  17),  ft  ia  aeen  that  for  L/A  *  0.4,  the  unmatched  broadeide 
gap  impedance  Zg(0,  0)  ia  capacitive]  for  L/A  -  0. 6,  it  ia  inductive)  while  for 
L/A  •  0. 5,  Z  (0,  0)  pasaes  through  resonance. 

8. 1. 3  FREQUENCY  AND  SCAN  DEPENDENCE 

Knowledge  of  frequency  and  scan-dependent  gap  impedance  Zg(  *,  k#Qj  f)  is 
needed  to  design  a  mulching  network.  Figure  18  illustrates  the  frequency  and  scan 
dependence  or  the  E-  and  H- plane  gap  impedance.  The  array  parameters  at  cen¬ 
ter  frequency  ffi  areb/Ac  ■  0.8,  d/Aft  ■  0.6,  s/Ac  •  0.28,  L/Ac  *  0.8,  h/Ac  ■  0.01, 
and  w/Ac  «  0.08  with  ka  «  120.  The  active  gap-impedance  curves  Z  (?,  k>Q)  are 
presented  for  the  three  frequencies  f  ■  0. 9fc,  f  »  f 0  and  f  •  1. 1  fQ.  In  Figure  16, 
the  gap  impedances  (normalised  to  ZQ  >  100 X?)  are  plotted  veraua  elevation  acan 
angle  0Q  and  discrete  values  of  asimuthal  acan  angles  *0*  where  coa  •o  "  ks0/k 

and  sin  0Q  sin  d0  «  v/ka,  »  -  0,  l . N  -  1.  This  enables  one  to  compare  the 

scan  dependence  of  the  dipole  gap-impedance  in  a  cylindrical  array  with  an  equiva¬ 
lent  infinite  planar  array.  Although  Zg  depends  on  both  the  scan  angle  and  fre¬ 
quency  f.  it  is  aeen  that  its  variation  with  dQ  is  much  greater  than  with  f  within  a 
±10  percent  frequency  band.  In  a  cylindrical  phased  array,  the  large  acan  angle 
dQ  ia  associated  with  the  edge  elements  of  the  excited  arc.  For  broadside- 
matched  elements  and  a  broadside  scanned  array,  these  elements  will  experience 
a  significant  impedance  mismatch.  However,  the  associated  reflection  power  loss 
ia  usually  reduced  by  an  amplitude  taper  of  the  excitation.  The  variation  of  gap 
impedance  versus  ($Q,  dg)  and  versus  frequency  ia  similar  (as  one  would  expect) 
to  that  of  a  corresponding  planar  array.  Figure  19  compares  gap  impedance 
Zg(*Q  *  n/2 ,  dg)  dependence  on  dg  at  f  »  fc,  for  ka  ■  30  and  ka  *  120  with  those 
of  the  equivalent  planar  array.  The  Smith  chart  normalisation  is  in  this  case 
ZQ  »  500 S3.  A  region  of  overlap  (dg  ■  CT  to  5<F )  may  be  noticed.  For  dg  >  8<F 
curves  for  ka  •  30  first  and  subsequently  those  for  ka  *  120  gradually  depart  from 
the  planar  ones,  due  to  the  curvature  effect.  Ik  was  verified  that  in  the  B-plane, 
the  values  of  Zg(  0Q,  d  >  ■  0* )  are  essentially  equal  to  those  of  the  active  impedance 
of  a  corresponding  planar  array  down  to  ka  *  30  and  are  therefore  not  shown. 

Figure  20  illustrates  the  frequency  dependence  of  the  H- plane  element  pattern. 
The  array  was  broadside -matched  at  center  frequency  fc.  A  simple,  realistic 
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Fleur*  17.  Gap  Impedance  Z«<0,  0)  Varsua 
Ground  Plan*  Spaefog  (b/X  »  u.8,  d/X  *  0.7, 
ka  ■  ISO,  Zo  *  lOOfi).  Parameter:  dipola 
length  L/X  ■  0.4,  0.5,  0.6 


1.0 


•1.0 


Figure  18.  H-Plane  Gap -Impedance  Versus 
Azimuth  Scan  Angle  and  E-Plane  Gap  Imped¬ 
ance  Varans  Elevation  Scan  Angle  Over  *10 
Percent  Frequency  Band  (b/Xc  »  0.5,  d/Xc 
•  0.8,  ka  -  ISO,  Z0  -  1001?) 
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Figure  19.  CWWiriloA  of  Cylindrical  and 
Planar  PMi  Amy  Gap-Impedance  for  ft* 
Plana  Scan  (b/rc  »  0.57d/Ac  *  0.  6,  20 
•  5000) 


Figure  90,  H-Pl*o*  Voltage  Fitment  Cain 
Patient  ft»/Ae  *  M»  d/*P  *0,1,  to  *  1?|), 
Parameter:  napwey  f  •  p.  9  fCf  rc,  1.  i  fc 
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matching  network  that  consists  of  a  A_/4  transformer  is  also  shown  in  Figure  20. 
The  transformer  is  placed  a  distance  away  from  the  dipole  so  that  the  trans¬ 
formed  gap  impedance  at  terminals  (2  -  2')  is  real.  Since  for  the  given  geometry 
Zg(0,  0)  *  (137  -  j38)42  and  Zj  =  50 Q,  the  value  of  Zg  was  chosen  to  be  Zg  =  13712. 
As  a  result,  4g  *  0.38^  and  Z^  =  9512.  With  these  (fixed)  matching  network  para¬ 
meters,  the  frequency  was  varied  over  a  ±10  percent  band.  It  was  observed  that 
the  voltage  element  gain  pattern  does  not  change  significantly.  The  pattern  shape 
depends  primarily  on  b /A  as  already  discussed,  while  broadside  gain  decreased 
0.2  dB  at  the  end  of  the  band.  The  curves  in  Figure  20  can  also  be  obtained  (in  the 
lit  region)  approximately  from 

lF  Ij".  <£)|  =  ”\|{1  -  |  r  1  (  0)|*  )  cos  <*>  (37) 

where  maybe  computed  from  Eq.  (19)  with  the  help  of  Figure  18. 

An  additional  remark  is  warranted  here.  It  is  known  that  in  planar  dipole  ar¬ 
rays,  the  feed  structure  can  give  rise  to  E-plane  blind  spots.  These  are  not  pre¬ 
dicted  by  analyses  that  disregard  the  influence  of  the  feed  structure  on  radiation. 
By  inference,  such  blind  spots  are  also  expected  to  appear  in  cylindrical  dipole 

arrays.  Various  empirically  based  remedies  can  be  employed  and  have  been 
o 

found  effective. 


5.2  Triangular  Grid 

Two  isosceles  triangular  grid  geometries  were  considered.  Both  grids  are 
derived  from  a  rectangular  lattice  with  b Ik  =  0.  6.  d/X  =  0.7,  and  preserve  the 
unit  cell  area.  The  first  grid  (a)  is  obtained  by  displacing  alternate  rings  by  half 
the  circumferential  spacing,  the  second  (b)  by  shifting  alternate  columns  by  half 
the  axial  spacing.  The  two  grids  produce  quite  different  element  patterns. 

In  the  isosceles  triangular  lattice,  for  a  constant  0  cut  the  EGL  drop-off,  due 
to  an  (m,  n)  grating  lobe  arriving  at  the  unit  circle  should  appear  at 


sin 


<t>  =  -T-fl  l  “~T  J  1-(cos  0  +  n 

sot  sin#  I  ml  N  \  26/ 


(m  =  0,  ±1, . . . ),  (n  =  0,  ±1, . . . ) 


(38) 


where  (m  +  n)  is  even. 
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9.  Herper,  J.  C. ,  Esposito,  F.  J. ,  Rothenberg,  C.,  and  Hessel,  A.  (1977)  Sur¬ 
face  resonances  in  a  radome  covered  dipole  array,  IEEE  Antennas  Propag. 
Symp.  Dig. ,  Stanford  University,  Stanford,  Calif.,  198-201. 


Figure  21  shows  the  local  grating  lobe  diagram  for  the  grid  of  type  (a).  It  is 
seen  that,  in  this  case,  relevant  grating  lobes  are  (-1,  -1)  and  (0,  -2).  According 
to  Eq.  (38),  an  off  H-plane  grating  lobe  (-1,  -1)  should  cause  narrowing  of  the  ele¬ 
ment  pattern  with  decreasing  9 .  This  is  in  contrast  to  the  rectangular  grid  refer¬ 
ence  array  Case  where  the  EGL  effect  is  due  to  an  H-plane  (-1,  0)  grating  lobe. 

As  a  result,  the  rectangular  grid  element  pattern  widens  with  9  tending  toward  the 
axis  [see  Figure  9  and  Eq.  (36)].  Based  on  Eq.  (38),  one  expects  for  b/X  =  0.6. 
d/X  -  0.7  that  #gQL  =  7S*  for  9  -  Q(f  and  =  51°  for  ^  =  7<T .  The  smooth 

pattern  in  the  H-plane  and  narrowing  of  the  patterns  (with  decreasing  9)  as  well  as 
the  predicted  values  of  are  also  observed  in  Figure  21.  In  addition,  one 

finds  an  increased  sharpness  of  the  EGL  drop-off  and  an  increased  level  of  the 
associated  off-broadside  ripple  with  decreasing  6.  Such  behavior  is  not  unexpected 
because  of  the  increase  of  mutual  coupling  when  the  grating  lobe  propagation  direc¬ 
tion  approaches  the  principal  H-plane  with  decreasing  9 . 

As  expected,  the  E -plane  element  gain  performance  is  found  to  be  similar  in 
both  the  triangular  and  the  rectangular  lattice  arrays,  since,  in  both  cases,  the 
axial  grating  lobe  is  located  at  -X/d. 

The  alternative  grid  (b)  and  the  pertinent  grating  lobe  diagram  are  shown  in 
Figure  22.  The  9  =  90“  pattern  is  similar  to  that  for  the  rectangular  grid  refer¬ 
ence  array,  which  is  plausible  because  of  the  same  grating  lobe  location.  In  the 
range  9  -  65°  to  90” ,  the  influence  cf  the  grating  lobe  (-2,  0)  on  the  element  pat¬ 
tern  is  dominant.  According  to  Eq.  (38),  however,  at  9  =  60* ,  two  EGL  effects 
are  expected:  one  for  0EGL  s  670  ^ue  to  scanned  off  H-plane  grating  lobe  (-2, 
0);  the  other  at  0EGL  3  32°  due  to  the  off- principal -plane  grating  lobe  (-1,  -1). 

Both  are  barely  noticeable  in  Figure  22.  As  9  decreases,  the  former  grating  lobe 
moves  toward  larger  values  of  0,  while  the  latter  approaches  0  *  (f  .  For  9  *  S0°, 
the  dip  is  predicted  to  appear  at  0EGL  =  18°  and*  toT  ®  ”  40°*  at  0.7®,  which  is 
well  borne  out  in  Figure  22.  For  a  further  decrease  in  9,  in  addition  to  (-1,  -1), 
the  (1,  -1)  grating  lobe  also  moves  into  real  space,  both  at  0  ■  0° .  This  results  in 
a  sharp  E-plane  dip  shown  for  9  ■  3 CP  .  With  increasing  values  of  0(for  9  ■  30° ), 
the  (1,  -1)  grating  lobe  moves  out  of  real  space,  which  produces  the  peak  in  the 
element  pattern  gain  near  2  CP  in  Figure  22. 

The  E-plane  gain  roll-off  is  slower  than  in  the  rectangular  grid  reference  ar¬ 
ray  case,  because  the  E -plane  grating  lobe  is  further  removed  from  the  real  space 
edge. 

In  view  of  the  above  data,  the  choice  of  rectangular  versus  triangular  grid  is 
based  on  the  maximum  axial  scan-angle,  the  extent  of  the  excited  array  arc,  and 
the  desired  frequency  bandwidth.  As  the  numerical  results  show,  the  lattice  de¬ 
sign  of  large,  cylindrical,  uniformly  spaced  arrays  may  be  safely  based  on  the  lo- 
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Figure  21.  Voltage  Element  Gain  Pattern 
and  Corresponding  Local  Grating  Lobe  Dia¬ 
gram  [Triangular  Grid  (a)b/A  =  6.  6,  d/A 
=  0.7,  ka  *  120].  Parameter:  elevation 
angle  0  ■  30P ,  40* » . . . ,  9CP 


AZIMUTH  ANGLE,  t  (DEO) 


Figure  22.  Voltage  Element  Gain  Pattern 
and  Corresponding  Local  Grating  Lobe  Dia¬ 
gram  [Triangular  Grid  (b)b/A  =  1.2,  i/A 
*  0.35,  ka  *  120].  Parameter:  elevation 
angle  0  •  30*,  40*,...,  90* 
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cal  grating  lobe  diagram.  (This,  of  course,  excludes  any  blind  spots  due  to  feed 
effects. ) 

5.3  Comparison  With  Other  Results 

The  few  relevant  references  found  in  the  literature  are  as  follows: 

1.  An  asymptotic  Green's  function  in  the  presence  of  a  conducting  cylinder 
is  employed5  to  analyze  circular  arrays  of  axial  dipoles  in  a  rectangular  lattice. 
Comparison  of  numerical  results5  with  those  based  on  modal  analysis  leads  to  the 
conclusion  that,  although  qualitative  element  pattern  features  in  the  forward  region 
such  as  the  EGL  drop-off  and  the  off -broad  side  ripple  are  exhibited,  the  details 
differ  from  those  obtained  from  modal  analysis  in  the  following:  (1)  The  9  »  7  ST 
near  0*0°  average  level  is  higher  than  that  in  the  principal  H -plane,  which  is  at 
variance  with  the  present  results  and  with  those  for  planar  arrays;  (2)  The  E-plane 
element  gain  loss  at  0  *  6(F  (off  axis)  is  about  25  percent,  as  compared  with  the 
32  percent  roll-off  predicted  by  modal  analysis;  (3)  No  shadow  region  results  are 
presented;  (4)  The  matching  scheme  and  the  pattern  normalization  are  unclear  and 
no  absolute  element  gain  information  is  presented. 

fl 

2.  Another  result  is  found.  Here,  the  analysis  was  carried  out  by  an  inter¬ 
esting  boundary  layer  method,  using  the  parabolic  equation  approximation  applied 
to  uniformly  spaced  arrays  on  convex  cylindrical  surfaces  of  arbitrary  but  slowly 
varying  curvature.  The  results  relevant  to  this  report  are  contained  in  Figure  3 
of  Indenbom  and  Filippov6  for  b/A  *  0. 6  and  d/k  =  0. 54.  For  convenience,  this 
figure  is  reproduced  in  Figure  23  here  along  with  results  for  6  *  60°  and  9(f  ob¬ 
tained  by  modal  analysis.  The  normalized  voltage  element  gain  pattern  in  the  6 

=  90°  cut  compares  well  with  that  based  on  modal  analysis.  However,  the  pattern 
for  0  =  60°  does  not  match  our  results,  both  in  the  0*0“  gain  level  and  in  tpe 
element  pattern  shape,  ft  exhibits  two  types  of  ripples  absent  in  our  results  and 
not  accompanied  by  the  usual  EGL  dips.  From  the  grating  lobe  structure  for  the 
geometry  in  question,  EGL  effects  are  predicted  in  this  case,  and,  therefore,  the 
pattern  is  expected  to  be  smooth, 

3.  We  have  also  compared  our  results  with  results  for  cylindrical  arrays  of 
aperture  elements.  The  mutual  coupling  effects,  such  as  the  endfire  grating  lobe 
condition  gain  drop-off  and  the  off-broadside  ripple  in  dipole  arrays  above  cylin¬ 
drical  ground  are  found  to  be  considerably  weaker  than  in  similar  arrays  of  aper- 
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ture  elements,  in  various  analyses. 1* 4* 10,  llf  13  This  is  expected,  in  view  of  the 
appearance  of  dipole  "images"  with  phase  reversed  currents  that  tend  to  reduce 
the  mutual  coupling  in  dipole  arrays.  By  contrast,  in  aperture  arrays,  the  "im¬ 
ages"  reinforce  the  excitation  of  mutual  coupling  effects.  It  is  reasonable  to  com¬ 
pare  the  element  patterns  for  arrays  of  axial  dipoles  with  those  of  circumferen¬ 
tially  polarized  rectangular  waveguide  elements,  since,  in  both  cases,  the  circum¬ 
ferential  plane  corresponds  to  the  direction  of  the  strongest  mutual  coupling.  It 
can  be  seen  from  Sureau  and  Hessel*0  and  Figure  24  that,  in  the  case  of  rectangu¬ 
lar  aperture  elements,  the  EGL  effect  appears  as  a  pronounced  dip  and  not  merely 
as  pscttern  drop-off.  Also,  the  off-broadside  ripple  level  is  appreciably  higher 
than  in  dipole  arrays. 


6.  CONCLUSION 

The  numerical  results  for  the  element  pattern  of  an  axial  dipole  element  in  a 
cylindrical  phased  array  lead  to  the  following  conclusions: 

1.  The  grating  lobe  effects  on  the  element  gain  and  phase  pattern  in  uniformly 
spaced,  large  radius,  circular  cylindrical  arrays  of  axial  dipoles  are  reliably 
predicted  by  locally  planar -grating  lobe  diagrams  for  either  rectangular  or  tri¬ 
angular  lattices. 

2.  The  creeping  waves  interfere  with  the  direct  element  radiation  in  the  for¬ 
ward  region  and  produce  a  far  field  in  the  shadow  region  with  an  angular  attenua- 

1  /  3 

tion  constant  proportional  to  (ka  sin  $)  ,  The  influence  of  both  effects  on  the 

element  pattern  (ripples  in  the  forward  region  and  finite  values  of  the  pattern  in 
the  shadow  region)  must  be  known  to  predict  the  sidelobe  performance  of  cylindri¬ 
cal  phased  arrays  accurately. 

3.  For  a  given  lattice  and  frequency,  the  element  gain  pattern  is  insensitive 


10.  Sureau,  J.  C. ,  and  Hessel,  A.  (1972)  Realized  gain  function  for  a  cylindrical 

array  of  open-ended  waveguides,  in  Phased  Array  Antennas,  A.  A.  Oliner 
and  G.H.  Knittel,  Eds.,  Artech  House,  Dedham,  Mass. 

11.  Borgiotti,  G.  V. ,  and  Balzano,  Q.  (1970)  Mutual  coupling  analysis  of  a  con¬ 

formal  array  of  elements  on  a  cylindrical  surface,  IEEE  Trans.  Antennas 
Propag. ,  AP-18(No.  1). 

12.  Borgiotti,  G.  V. ,  and  Balzano,  Q.  (1972)  Analysis  and  element  pattern  design 

of  periodic  arrays  of  circular  apertures  on  conducting  cylinders,  IEEE 
Trans.  Antennas  Propag. »  AP-20(No.  5). 

13.  Borgiotti,  G.  V. ,  and  Balzano,  Q.  (1972)  Conformal  arrays  on  surfaces  with 

rotational  symmetry,  in  Phased  Array  Antennas.  A.  A.  Oliner  and  G.  H. 
Knittel,  Eds. ,  Artech  House,  Dedham,  Mass. 

14.  Collin,  R.E.  (1960)  Field  Theory  of  Guided  Waves,  McGraw-Hill. 
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Figure  23.  Comparison  of  Voltage 
Element  Gain  Patterns  Between 
Figure  3®  and  Results  Based  on 
Modal  Analysis  (b/X  *  0, 6,  d/X 

*  0.54,  WX  •  Q.  5.  s/X  ■  0.25, 
w/X  •  0.05,  h/X  *  0.01,  ka 

*  125.7).  Parameter!  elevation 
angle  $  •  80P»  90° 


0  20  40  60  80  100  120 

AZIMUTH  ANQLE,  0  (DEO) 


Figure  24.  Comparison  of  Voltage  Element  Gain 
Patterns  Between  Dipole  and  Rectangular  Aper¬ 
ture  Radiators,  ^  (b/X  *  0,6,  d/X  *  0.8,  ka  *  105) 
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to  the  element  dimensions  and  ground  plane  spacing,  provided  in  each  case  that 
the  element  is  appropriately  "broadside"  matched. 

4.  The  unmatched  gap  impedance  exhibits  significant  dependence  on  dipole 
length  and  its  distance  to  ground.  Thus,  the  basic  limitations  of  the  dipole  ele¬ 
ment  rests  in  its  scan  and  frequency-dependent  impedance  properties  and  those  of 
the  matching  scheme. 

5.  The  gap  impedance  is  similar  to  that  of  a  corresponding  planar  array. 
Therefore,  planar  array  simulator  techniques  for  designing  and  matching  radiat¬ 
ing  elements  can  be  used. 

6.  The  dii>ole  phase  center  is  located  near  the  element  and  not  on  the  ground¬ 
ing  surface. 

7.  Mutual  coupling  effects,  except  for  possible  E-plane  feed  interactions  pro¬ 
ducing  blind  spots,  are  less  pronounced  than  in  similar  aperture  arrays. 

8.  Polarization  purity,  wide  bandwidth,  lower  mutual  coupling  effects,  ease 
of  fabrication,  and  low  cost  make  the  axial  dipole  a  suitable  element  for  cylindri¬ 
cal  phased  arrays. 

9.  The  theoretical  results  compare  well  with  experiment  (see  Part  2  of  this 
report). 
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Appendix  A 


Floquet's  Theorem  for.  Isosceles  Triangular 

Array  Lattice 


With  reference  to  Figure  4*  let  die  usual  Floquet  expansion  for  any  field  com¬ 
ponent  in  the  unit  cell  of  sise  (2rr/N  x  id  )  he 

(Ala) 

where  amjft  are  the  amplitude  coefficients  of  the  (m,  n)th  Floquet  mode,  with 

v  =  v+  mN  =  v  +  ~  m  (m  ■  o,  ±1, . . . )  (Alb) 

m  Q 


and 


*4ns *«o  +  d  n 


(n  *0.  ti»  i  •  •  )• 


For  a  triangular  lattice.  Floquet's  theorem  imposes  an  additional  constraint 
on  Eq.  (Ala)  so  that  a  displacement  by  (ir/N)  in  ^  followed  by  that  of  d  in  z  multi¬ 
plies  the  field  by  the  phase  constant  exp(-jd).  that  is. 


¥{*  +  £.  *+d)  =a-J*¥(<M). 


(A2) 
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It  is  seen  from  Fleur*  A1  that 


(AS) 


Flfur*  Al.  Unit  Cell  Cross  Section  at 
P  ■  PQ  of  a  Circular  Cylindrical  Array 
With  Triangular  Lattice 


Substitutlne  Eq*k  (Ala)  and  (AS)  into  Sq.  (Al),  on*  finds 

5  •**'">  #  4  k*"d> 

_  ^  Ww 

m,n*  '* 

=  +  ^  .  (A4) 

Using  Eqs.  (Alb)  and  (Ale)  and  the  orthogonality  of  the  Floquet  basis  in  Eq.  (Ala), 


a-iKkf*  ^(imrenw) 


(A5) 


(AS) 


whereupon,  it  follows  that 
a|nn(*-Hn.*n)ir_1,  =  0 
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Thus.  «mn  »  0  ttnlaaa  (m  +  n)  la  an  awn  lntttar.  Consequently,  we  conclude 
that  la  an  iaoacalaa  triangular  lattice*  Floquat  modes  with  (m  +  a)  odd  are  not  ex¬ 
cited. 
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Appendix  B 

Single  Element  Excitation  in  a  Triangular  Lattice 


A  single  (s  •  0.  t  ■  0)  element  excitation  in  the  triangular  isosceles  lattice  of 
Figure  4  is  expressed  in  terms  of  a a  excitation  vector  a  •  f  (a  *  0,  1. , . . , 
N  -  U  t  -  0,±1. . . .)  where 


d 

47TN 


N-1 
V 

V*~M 


i  a 


2 £» 
a 


ir/d 

/ 


JkNdt 


«*»o 


-sr/d 


©1) 


Each  exponential  term  in  Eq,  ©1)  represents  a  unity  -amplitude  progressively* 
phased  array  excitation  2*WN  in  the  angular  direction  and  k>Qd  +  W¥f  N  in  the 
skew  direction.  In  this  Appendix,  we  verify  the  identity  Eq.  ©1).  To  this  end, 
one  observe 8  that 


ir/d 


-ir/d 


(t«0.±l. ...)  . 


©2a) 


1 

ir 


“In 


2  [•**  " 


I 


I  skl7T(2t+t) 

*  n*  Z*™  9  *  **  **"  *  **"°* 11  *****  ’ 

The  product  of  the  loft-hand  aidoa  of  Eqs*  (B2a)  and  (B2b)  la 
$to  fa  *  ®to  ^«.tn* 


(B2b) 


(BS) 


Since  a  *  0.  1. . . . .  N  -  1,  it  follows  that,  in  Eq.  (B3),  only  the  1  •  0  term  contri 
butes.  Th.r.for«.  4(o  +,  ,1N  ■  l,0  ^Q.E.D. 
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Appendix  C 
Covergence  Acceleration 


As  already  mentioned  in  Section  4,  the  double  sum  in  Eq.  (16b)  may  be  evalu¬ 
ated  as  £  gQ  (£  fmn).  In  this  Appendix*  we  describe  a  convergence  acceleration 

of  the  series  £  r  To  this  end*  we  write 
m  mn 


an~a~  2  fmn 


where  f|nn  is  given  by  Eq.  (16b)  as 


(Cla) 


*in*(i»miO 

C-mX)2 


(m  *  0,  ±1* » . . ; 


Z|rm(Kl»;  •*  p)  (*nPo'  n  a  o,  ±1.  ...  ) 


(Clb) 


with  x  -  w/2p  .  For  convenience*  let  <7  =  a  +  a',  where  o  ■  f  _  and 
o  o  o  on 


(C2) 


The  prime  indicates  the  exclusion  of  the  m  =  0  term  in  the  infinite  sum.  The 
method  to  be  described  is  carried  out  for  both  the  rectangular  and  the  triangular 
lattice  geometries. 
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Cl.  RECTANGULAR  LATTICE 


To  facilitate  convergence  acceleration  of  o',  it  is  desirable  to  introduce  the 

(2) 

first-order  large  index  asymptotic  expression  for  Z^H ym,  that  is. 


(2)  1 
Zi>m  (*n  ?  •*  P)  Hvm  ^nPo)  “  *  J  Yi»m  <*nPo>  ~  J  ^ 


7TlmlN 


(C3) 


In  vie-w  of  Eq.  (C3),  one  can  rewrite  Eq.  (C2)  in  the  following  form: 


Fmn  +  ® 


m»-* 


(C4a) 


where,  since  u  *  u  +  mN, 
m 


•in2(vmx)  i 

F- - —  I V.  K  ■■  »•  Pa)  H®  |*„p0)  -  i 


'm 


,  *  -  s  i/*+2i/rrtN  ; 

°  ‘  n ImlN  ‘  1  jrlmPN5' 1  * 


and 


1  y  **2iVm*\ 
1  ttN3x2  mT-»  Iml3 


<C4c) 


o 

One  observes  that  the  series  in  Eq.  (C4b)  now  converges  as  1/m  .  The  stum  S  ih 
Eq.  (C4c)  may  be  evaluated  with  the  aid  of  Collin, 14  p.  S79.  The  result  is  the 
rapidly  convergent  series: 


S*j  1  IA(1  -  cos(2Px) )  - B] 


(C5a) 


where 


A  =  I  1 .20208690 

m»va.3..  .,IH 


(CSb) 


and 


B  =  (Nx)2  cos  (2i/x)  [2  In  (2Nx)  -  3  - 


(Nx>2 

18 


■  M« 


]  ‘ 


(C5c) 
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C2.  TRIANGULAR  LATTICE 

In  view  of  Eq.  (A6),  we  distinguish  two  cases: 
For  n  odd 


1.3,... 

a-  £  Fmn  +  Sodd  (n  ■  ±1.  ±3. . . . ). 

n=...,-3.-i. 


The  series 


14 

may  be  cast  into  a  rapidly  convergent  form  (Collin,  p.  580) 


sodd  =  J  lAodd<1-co«(2vx)-Bodd] 


where 


-  1 

Aodd=  2  —3  *1.05170979 

m=i.3....  m 


and 


a  (Nx)2 

Bodd  =  (Nx)2  cos  (2 ux)  (In(Nx)  -  +  -rr~  -  •••  ]  . 


2  18 


When  n  is  even 

2.4. 

1 

m  =....-4.-2, 


0  =  ao  +  ^  Fmn  +  S, 


where 


Se  =  S  '  Sodd* 


(C6) 


(C7a) 


(C7b) 


(C7c) 


(C7d) 


(C8a) 


(C8b) 
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Appendix  D 
Matching  Networks 

The  load  impedance  is  connected  to  the  output  terminals  of  a  two-port  match¬ 
ing  network  characterized  by  its  scattering  matrix. 

In  Section  Dl.  we  first  derive  expressions  for  the  reflection  coefficient  at  the 
input  to  the  matching  network  and  for  total  voltage  at  the  load. 

Similarly,  in  Section  D2,  we  derive  expressions  for  the  reflection  coefficient 
at  the  input  to  the  matching  network  and  for  voltage  at  the  load  (dipole  gap-imped¬ 
ance)  which  is  now  connected  to  output  terminals  of  a  matching  network  via  a 
transmission  line. 

In  Sections  D3  to  D5.  we  determine  the  elements  of  a  scattering  matrix  for 
three  specific  simple  matching  networks:  series,  parallel,  and  quarter-wave 
transformer. 


Dl.  LOAD  IMPEDANCE  CONNECTED  TO  OUTPUT  TERMINALS  OF 
TWO-PORT  MATCHING  NETWORKS 


The  S-parameters  of  a  linear,  passive  two-port  network  can  be  defined  with 
the  help  of  Figure  Dl.  As  shown  in  the  figure,  vj”  (j  ■  1,  2)  are  the  traveling 
waves  incident  on  the  junction,  and  Vj  are  the  reflected  traveling  waves  from  the 
junction.  By  linearity,  the  incident  and  reflected  waves  are  related  by 


S..  V7 
il  • 


for  i 


1.  2 


(Dla) 
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Figure  Dl.  S -Parameter  Representation  of  a  Two-Port 
Network 


or 

VI  ■  Sll  V1  +  S12  V2  ®»*>> 

V2  ’  Sji  Vj  +  S22  V2  .  (Die) 

From  Eq.  (Dlb),  it  is  seen  that  is  the  reflection  coefficient  in  the  input  guide  1» 
with  guide  2  terminated  in  a  matched  load.  Also*.  S21  is  the  transmission  coeffi¬ 
cient  into  guide  2  from  guide  1  under  the  same  condition.  Similar  statements  ap¬ 
ply  to  the  parameters  S22  and  S^* 

If  guide  2,  with  its  characteristic  impedance  Z2.  is  terminated  by  an  imped¬ 
ance  Z22-  at  the  terminal  plane  (2-2')  (see  Figure  D2),  then  Vj  may  be  regarded 
as  the  incident  wave  on  Z22,,  and  V2  is  the  wave  reflected  from  Z22  The  ratio 
must  be  equal  to  the  reflection  coefficient  of  the  load;  hence. 


Z22-'  -  ^2 
Z22'+  Z2 


(t>2) 


Figure  D2.  Two-Port  Matching  Network  Connecting 
Load  Impedance  Z22"  With  Transmission  Line  1 
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Using  Eqs.  (Dl)  and  <D2),  one  can  now  determine  the  input  reflection  coefficient 
/}  in  terms  of  S  parameters  and  7^: 


v;  ■  su  vi  *  su  vi  ■  su  r2  va 

(D3a) 

si2vt"v2-saarjv2* 

(DSb) 

so  that 

r  =  VI  =  r  .  S12  sai  r2  M  Su  -  dr2 

1  11  ^22  r2  “  1  1  *  ^2  r2 

(D3c) 

with 

4  *  S11  S22  *  S12  S21* 

(D3d) 

Consequently,  the  input  impedance  is 

Zin  *  Zi  1  +  ^  • 
in  1  i-Fl 

(D4) 

One  can  also  find  total  voltage  Vg  at  the  load  in  terms  of  the  S  parameters  and 
the  incident  voltage  V+: 

In  view  of  Eq.  (D2),  one  has 

V2_  •S21  V1  +S22  r2V2 

(D5a) 

or 

S21  4- 

(D5b) 

*  ■  1  -  Sj2  r2  V 

Using 

v"  =  r  v+ 
vi  i  vi 

(D6a) 
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-aliiii.:.;.. . • .  1.1 


one  finds  from  Eq.  (Dlb) 


*  Vi  -  Su  v}  n  Vj  -  Si!  vt  n  -  Sn  + 

Vo  - - s - * - s - *  T - V1 

2  S12  S12  S12  1 


(D6b) 


or,  upon  substitution  of  (/"j  -  Sjj),  from  Eq.  (D3c) 


v+  a  _  sji  r,2«  v+ 

2  ®22  * 


(t)6c) 


Using  Eqs.  (D5b)  and  (D6c),  the  tdtsl  voltage  at  (2  -2')  is 

+  -  s21  (i  +  r2>  + 

v2  a  v2  +  v2  *  i  -  S22  r2  V1* 

When  the  input  line  is  matched  as  shown  ih  Figure  t>3,  one  readily  sees  that 


V2aV2+V2isV2(1+S22) 


(D8a) 


+  V2 
V2  =  1  +  S„ 


(D8b) 


and,  consequently*  since  =  0, 


+  s12 

r  1  - 


VlaVl=S12  V2aTTs^  V2- 


(bdc) 


The  output  impedance  is  then 


Zout  =  Z2 


1  +  S22 

1  -  Soo 
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Figure  D3.  Two-Port  Network  Connecting  Matched- 
Transmission  Line  1  With  Transmission  Line  2 


D2  LOAD  (DIPOLE  GAP  IMPEDANCE)  CONNECTED  TO  OUTPUT  TERMINALS 
OF  TWO-PORT  NETWORK  VIA  TRANSMISSION  LINE 

Transmission  line  of  length  ig  is  terminated  with  dipole  gap  (load)  impedance 
Zg  and  connected  to  output  terminals  of  matching  network  as  shown  in  Figure  D4. 
The  output  transmission  line  2  has  characteristic  impedance  Zg  and  propagation 
constant  02* 


In  this  case,  from  Eq.  (D3c),  the  reflection  coefficient  at  the  input  to  the 
matching  network  is 


SS 


(DlOa) 


S11  ^ 

1  -  S22  r2 


where  now 


r2  -  rg 


(DiOb) 


Z«  -  Z2 


'■VS 


(DIOc) 


id  total  voltage  at  z  -  lj,  that  is.  Vg  *  V<z  *  Aj).  interns  of  incidence  volt' 
j  at  the  input  port  (1  -  1")  of  the  matching  network,  it  ia  convenient  to  pro- 


To  find  total  vol' 
age  V  j  at  the  in] 
ceed  as  follows: 


V2  (z  =  0+)  «  V2  (z  «  0+>  +  V2  <»  •  0+)  »  Vj  (a  -  0+)(l  +  r2)  (Dlla) 


where 


2  Vj  (z  «  0+>  * 


miib) 


and  rg  is  given  by  Eq.  (DIOc),  Using  relations 
V2(z)  =  V2  (z  *  0+)  e"^2*  +  V“  (z  =■  0*>  e1^2* 


*  vj  <z  =  o+)  e‘^a*  a  +  r2  e^2*) 


(D12a) 


and  Eq.  (Dllb),  total  voltage  at  the  load  is 

v2  U  3  V  =  Vg  *  V2  <z  a  0+)  e'^2i2  <1  +  rg). 

Substitution  of  Eq.  (D12a)  into  Eq,  (D12b)  yields 


V  -  V(,.lVA‘!  • 


(Diab) 


(DIS) 
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Furthermore,  substituting  Eq.  (D7)  into  Eq.  (D13),  we  finally  obtain 


S21(l+/g )o.1^l0  v+  . 

i  -  s22  r2 


(D14) 


D3.  S-PARAMETERS  OF  THE  SERIES  MATCHING  NETWORK 

The  scattering-matrix  elements  will  be  evaluated  for  the  series -matching  net- 
work  connecting  two  lines  with  characteristic  impedances  Zj  and  Z2  as  shown  in 
Figure  D5.  This  matching  network,  which  consists  of  transformer  and  series  re¬ 
actance  in  cascade  connection,  .is  employed  in  Section  D5  for  matching  the  dipole 
gap  impedance  Z^(v,  kzQ)  to  TEM  feed-transmission  line. 


The  values  of  the  transformer  ratio  Dq  and  the  reactance  X22*(v,  kz0)  are 
chosen  so  that  active  reflection  coefficient  r^(v.  kZQ>  at  the  input  side  of  the 
matching  network  is  zero.  To  determine  Dq  and  kzQ),  we  first  write  the 

expression  for  impedance  at  the  output  terminals  (2  -  2  ")  of  the  matching  network. 
In  reference  to  Figure  D4,  the  impedance  looking  into  transmission  line  2  at  z  -  0 
is 


V  *  ^zO^  3  ^z(P  +  3  ^22*'^*  ^*0^ 

Zg  (v,  kzp)  cos  02*2  +  i  z2  sin  &2l2 
~  Z2  Z2  cos  P2l2  +  jZgU»  kz0)  sin^242 

where  the  gap  impedance 


Zg(v,  kZQ)  *  Rg(v,  kjg)  +  j  I^q). 


(D15a) 


(D15b) 
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Thus. 


no  a  V“ 


1*22  *«<)> 


(D16) 


With  elements  of  the  series  matching  network  so  defined,  the  corresponding 
S-parameters  can  be  found  as  follows: 

If  the  output  line  is  matched  as  shown  in  Figure  D6,  we  have 


Zg -_jX2a^»,  k,Ql 

V£  I  b  Sin  ■  Z1  ^  _ _ _ 

vt  I  "  zin  +  Z1  ‘  Z2“3x22-^»^q),  „ 

tr.A  *  r 


Z  2  ~  Zja^IVi  kto^ 
*  e 

Z2  +  z22'(**  ^sO^ 


where  Z^(y,  kz0)  is  given  by  Eq.  (D15a). 


<P17> 


Figure  D6.  Network  Illustrating  Evaluation  of  Sjj 
and  Sa  i  of  Series  Matching  Network 


With  the  input  line  matched  (see  Figure  D7),  one  has 


Zout  -  Zg 


Zout  +  Z2 


_  Z22  "  (y»  kzQ)  -  Za 
Z2  +  Z22'^v*  ^*0^ 


2 

2inp  -  jX22"(y.  kz0)  -  Z2 
Zln0  "  &22'^v>  +  Z2 


(D18) 
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1 

i  1:n*  -JXtr 

1 

2 

1 

1 

1 

1 

v}2*’>a 

T 

i 

1 

_  J 

2* 

Figure  D7.  Network  Illustrating  Evaluation  of  Sgg 
and  Sj2  of  Sariaa  Matching  Network 


To  find  S3  j,  we  again  consider  the  output  line  matched  (see  Figure  D6) 
the  input  line,  we  have 

vi '  vi +  vi  ■  vi(l  *  sn* 


and 


ll  *  Y1  <V1  “  V?  *  Y1  V1  <l  “  *li* 


Since 


h  Ylvl 


and 


■i  *  Ya  va 


we  see  that 


Y2  Vi  " 


YiVf 


n. 


(1  -  sn). 


We  now  obtain 


Yl  Z-2 

- (1  -  s  )  - - - - 

Y2n0  Z2  +  i^O  Zj  **  J  x22'(|,»  kz0) 


Zg  2  kept  Zi 

3  zl  z2+Z22'{Vt  k*0} 


.  On 

(D19a) 


(D19b) 


(D20a) 


(D20b) 


(D20c) 


(D21) 
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Dk  S-PARAMETEHS  OF  A  PARALLEL  MATCHING  NETWORK 

Her*,  we  evaluate  a  scattering  matrix  of  a  parallel  two  port  matching  network 
which  consists  of  the  transformer  and  susoeptance  in  cascade  connection*  as  shown 
in  Figure  D8.  . 

The  values  of  the  transformer  ratio  nQ  and  the  susoeptance  (a.  kgQ)  are 


Figure  D8.  Parallel  Matching  Network 


choaen  so  that  tha  act  Its  reflection  coefficient  ^  (»•  k#Q)^  tha  input  aida  of  tha 
matching  oat  work  ia  saro  for  tha  apaeifle  aat  of  r ,  k^  (aaa  Figura  D4).  Tha 
tranaformar  ratio  n^  and  sosceptanee  ®2  j  -  (*.  k#Q)  can  ba  datarminad  from  tha 
ralation  for  admittanca  at  tha  output  terminals  of  tha  matching  network  (2  -  2') 
looking  into  transmission  line  2,  which  is 

k*°’  ■  z  77T~i  •  0M'(  '•  k.o>  *  *11'  k.0>  <DJ«> 

k,o> 


whara  Z]2«  (a,  k#Q)  ia  glean by  Eq.  (D13a).  It  is  saan  that 

n°  h  • 


(D25) 


Next,  wa  determine  tha  S-paraxneters  of  tha  parallel  matching  network.  With 
tha  help  of  Figure  DS,  wa  write 


*11 


*1  -  *1  •  d  Oft  -  W> 

,  Y»  +  Yl.  Y1  ♦’4,T*  -  J8!*-''-  k,0>» 

*  A 


Vj-0 


kjO )  -  Yj 
Yjj  *(r,  +  Yj 


and 


<D26a> 


c  V3 
V1 


vj-o 


Yj  2Yj  l 

- (1  -  S,,)  »  —  - 

Ya°o  “  «o  Yia'“’-k.o>  +  Ya 


(D2flb) 


With  tha  aid  of  Figura  DIO,  following  a  similar  procadura,  one  determines 
S22.  that  is. 
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Y2  -  Y2f>(y,  kz0) 
^2  +  Y22'*^I/*  kzO^ 


(D37a) 


1  I  1:nG  -Tb22'  T  2 


Figure  D9.  Network  Illustrating  Evaluation  of 
Sjj  and  S21  of  the  Parallel  Matching  Network 


Figure  DIO.  Network  Illustrating  Evaluation  of 
S22  and  S12  of  the  Parallel  Matching  Network 


while  S12  can  be  determined  from  the  same  figure  following  the  procedure  indi¬ 
cated  in  the  previous  section  or  via  a  reciprocity  relation,  that  is, 

S12  '~S31  •  <DJ7b> 

M 
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DS.  S-PARAMETERS  OF  THE  QUARTER-WAVELENGTH  TRANSFORMER 
MATCHING  NETWORK 

A  matching  network*  shown  in  Figure  Dll,  consists  of  a  transmission  line 
with  characteristic  impedance  and  propagation  constant  The  transmission 
line  length  is  =  A^/4.  The  network  is  placed  a  distance  I2  away  from  the  dipole 
so  that  the  transformed  gap  impedance  Z22*  at  terminals  (2  -  2"*)  is  real,  that  is, 

**\Z22Av'  kzO;Vls0  (D28) 

where  Z^-  * V •  ^0^2^  is  8iven  by  Eq.  (D15a). 


Figure  Dll.  Quarter-Wave  Transformer  Matching  Network 


From  Eq.  (D28),  we  obtain  the  relation  for 
„  2Z2  Xg  ( v ,  kzg) 

*8  2^212  =  - T2 - 2 - ? 

R|(,.ki0)+Z2+X|(,,.ltzo> 

where  gap  impedance 

V-  V’V”-  V 

and  characteristic  impedance  of  transmission  line  2  is 

Z2  =  Rg  (v'  kz0)  * 


(D29a) 


(D29b) 


(D29c) 
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<D30a) 


The  characteristic  impedance  of  A^/4  transformer  is  then 

2t  *  '^ZiRetZ22'(l''  ks0*  *2)] 
and  the  propagation  constant 

(D30b) 

2 

In  the  remaining  part  of  this  section,  we  determine  the  scattering  parameters 
for  a  section  of  transmission  line  of  length  X^  and  characteristic  impedance  Zt 
shown  in  Figure  D12. 
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!  *»  >p» 
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L  _  _ 

J  * 
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Figure  D12,  Quarter-Wavelength  Transformer  Matching 
Network  Illustrating  Evaluation  of  S-Paraineters 


Scattering  parameters  attd  S21  are  found  with  the  help  of  Figiire  DiS  ttS 
follows: 

We  define  reflection  coefficients 


z2  -  Zt 

t(z  =  Xt)  »  =  z2  +  Zt 


(bain) 


arid 


r(z  =  o+)  •  r2  e-J23*'* . 

The  impedance  at  the  iriput  port  of  the  matching  network  is 

1  +  t(z  =  Q'1') 

Z(z  =  0)  =  Zt  j  .piz  =  0+)  * 


(baitj 
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1  r - 

”1  2 

v: 

!  0  1 

! 

!  0  1 

> 

•  1 

1  29 

2*0” 

2-04 

Z< 

>0 

Z*i,  2 

Figure  D13.  Network  Illustrating  Evaluation  of 
Si  i  and  S*i  of  Quarter-Wave  Transformer 
fetwork 


Consequently,  we  may  write 


V! 


'11 


Z(z  *  0)  -  Z\ 


(D33) 


Vj  1  Z(z  *  0)  +  Zj 

vj=0 


Substituting  Eq.  (D31a)  into  Eq.  (D31b),  then  Eq.  (03  lb)  into  Eq.  (D32),  and 
finally  Eq.  (D32)  into  Eq.  (D33),  we  obtain 


'11 


<Z2  -  ZtXZt  +  Z|)  +  (Zt  -  Z xHZ2  +  Zt) 

(Z  2  -  Zt)(Zt  -  Zj)  e‘*4*t*t  +  (zt  +  ZX)(Z2  +  Zt> 


To  find  S2  j,  which  is  defined  as 


V2 


21  v+ 


(D35a) 


we  first  write  the  following  relevant  expression  for  the  reflection  coefficient  at 
z  =  0” : 


T(z  =  0")  *  » 


Z(z  *  0)  -  Zj 


Z(z  -  0)  +  Z 


(D35b) 


1 
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For  0  <  *  it,  total  voltage  on  the  transmission  line  is 

V(s>  »  V*<0+>  e"^*  (l  +~PT" 

\  V+(0+)  / 


At  *  «  lt, 


vut)  •  v+(o+)  e*3*4*  a  +no+)  J2^) . 


Substituting  Eq.  (D31b)  into  Eq.  (D37),  we  have 
+  +  "i&*t 

v(it) « v+(o+)  e  1  *  a  +  r2) . 

Furthermore,  since 

V(0+)  »  V(O')  «  V(0) 
we  see  that 


V<0+>  *  V+(0+Kl  +T(0+)) 


and  consequently,  using  Eqs.  (D39a)  and  (03  lb)  from  Eq.  (D3fib) 


+  +  YM. 

V  <0+) 


V(0) 


i  +  r<o+)  i  +  r«  e’32^4* 


Here 


V(0)  *  V«f)  -  V+(0')  (1  +  r:) 


which,  after  substituting  into  Eq.  (D39c),  yields 

±ft.> 

WMt  * 


.  .  v+(o*)d  +  r i ) 

V+(0  )  *  — - 


i  +  r2e 


(D36) 


(037) 


(038) 


(D38a) 


(038b) 


(0380 


(D40a) 


(040b) 
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Substituting  Eq.  (D40b)  into  Eq.  (D38),  wo  obtain 


vut) 


V~V)  (1  ♦  fj)  e"^*4*  (l  +  Ta) 


i  +  r„  o 


•m*x 


(D41) 


Using 

V(lt)-V“(s-iJ')  (D42a) 

one  can  therefore  write 


V‘(*  .A? )  (1  +  H)  (1  +  To)  -tf.*. 

-  ^77,  ■  ■,  +v2  ,W.  • 


Finally,  substituting  Eq.  (D35b)  for  /  j  into  Eq.  (D42b),  we  obtain 


S21  *  ‘ 


2Zt  (i  +  r2> 


-ifi±A 


t*t 


zt  (l  +  r2  e“Ja^tAt)  +  Zjd  -  r2  e‘32^iAt) 


(D42b) 


(D43) 


where  is  given  by  Eq.  (D31a). 

With  the  help  of  Figure  D14,  following  the  same  procedure,  one  finds  S22  and 
S12.  However,  because  of  the  symmetry  of  the  matching  network,  one  sees  from 
Figures  D13  and  D14  that  expressions  for  can  be  deduced  from  Eq.  (D34)  by 
interchanging  subscripts  1  and  2.  Thus, 


VJ  (Zt  -  ZtHZt  +  Z2)  e  +  (Zt  -  Z2)(Z j  +  Zt) 

V+  (Zj  -  Zt)(Zt  -  Z2)  e  1  +  (Zt  +  Z2)(Zj  +  Zt) 


(D44) 


Using  a  reciprocity  relation,  one  can  also  write 

S12  ■  r1 321  «d«> 

Z2 

where  S2  j  is  given  by  Eq.  (D43). 
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Figure  D14,  Network  Illuet  rating  Evaluation  of 
Sag  wnd  Sig  of  Quarter-Wave  Tranaformer 
Matching  Network 
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Appendix  E 

Asymptotic  Expansions  of  2„j («n*xnp^ 


In  this  Appendix,  we  derive  asymptotic  expressions  of  zrm(*Q**  KnP0  *  and 
zPm(*n*‘  *n^0>  H^m  z*rm<*na*  *n^>)  ^  «ivan  by  -Bq‘  <7b)  ** 


•)  fa) 

Z|'m(,tna«,fnP0>  *  ^m^nV  "  „  (2)  Hym  (*nV  * 


Una) 


(El) 


For  simplicity  of  notation,  we  write 


*1*  *n* 


*2  5  *nP0 
1,3  "m 


(E2a) 

(E2b) 

(E2c) 


where 

K  3  >/k2  -  k  2  ,  L[k  ]<  0. 
n  zn  •  *m‘*nJ 

In  this  view,  Eq.  (El)  takes  the  form: 


(E2d) 
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JyiXi) 


(E3a) 


2y(Xj  ,  Xj)  *  JpfXg)  - 


■«) 


r(2) 


<*2) 


<*1> 


where 

(x)  «  J„  (x)  -  j  Y,  (x).  <®«») 

The  asymptotic  expansions  of  Eq.  (E3a)  will  be  derived  in  all  relevant  regions  with 
respect  to  orders  and  arguments  for  berth  real  and  Imaginary  argument  i. 


El.  REAL  ARGUMENTS 

El.  1  Debye's  Asymptotic  Expansion  for  Large  Orders 

If  x  is  fixed  and  positive  and  v  is  large  and  positive,  then  from  Eqs.  (9. 3. 7), 
(9.3.8).  and  (9.S.9)7 

J„  (x)  -  -T-fcr  s>  <E4a) 

/2irv  tanho  J 


V,  (x)  -  -  ./■  !  .  ««b> 

Vxrtanha  7 


(Sic) 


(S4d) 

00  u v  (cotha) 

ii  1  1  +  t  -  (B4e) 

j  k«l  *k 

(cotha) 

3  *  1  +  I  (-l)k  T  .  (E4f) 

y  k»i  v11 


where 


P  *  a  -  tanha 

1  1  +  M  -  (f 


a  «  — £n 
2 


) 


1  -  ^1  -  (f)8 
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In  Eqa.  (E4e)  and  (B4f) 


«0  (t>  “  1 

tt  (t) « 3Li  ,*3 
1  24 

,tt.  »lta  -  4W  t4  ♦  »M  t» 
^  1152 


<E5a> 

(E5b) 

(ESc) 


u3  (t) 


u4  (t)  - 


30375  t3  -  369503  tS  +  748765  t7  -  428425  t9 
414720 


4465125  t4  -  04121676  t*  +  340022430  t8  -  446185740  tl0+ 185010725  t 


(E5d) 
12 


308  13120 


<E5a) 


or  In  ganaral 


t 

(1  - 18> ^  (I)  J  i\  -  5ta)  ^  tt)  dt.  (k  •  0.  !*...). 

2  8  n 


<ESt> 


In  terma  of  Beaael'a  function  of  fir  at  kind  J*  (x)  and  Neumann  functiona. 
Yv  (x),  Zv(xy  Xj)  can  be  written  in  the  form 


Zv(xj.  Xj)  «  Jp  (xj)  - 


1  rji<x1)Jp(x2)  +  Jp(xj)  Yp(xj)  Yp(x2) 


F(x,) 


[ 


Ya  (Xj) 


J»(X|)  Jyixg)  YpUj)  -  Jp  (X|)  Yp  (xj)  | 
Ya  (Xj) 


(E6a) 


where 


F  (x.)  -  1  +  _ 

X  y4 


j*  (*,> 


(E6b) 


y;  <*i} 
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1 


Substitution  of  Eq.  (E4)  into  £q.  (E«)  yislda 


**'  ’  /l„.«nko,  nj  F^.l.JL 

+  *yi«y»  -  A J>  +  ,-»<♦*  -  *»)  .  *y1*i»  .-»<«*» 


(E?a> 


where  asymptotic  expansion  of  F(Xj)  is 


V  ■'*(%) 


,*4»0l  . 


(E7b) 


In  Eq.  (El).  subscripts  1  and  8  rafar  to  arguments  Xj  and  Xj.  respectively* 
Similarly,  on*  may  writs 


from  where  using  Eq.  (E4),  the  rsspsetivs  asymptotic  expansion  is 
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1 


Z¥{xvx2)  H  ^Nxj) 


»*tanhas 


a 


+  i#ja  *ya 


1  /^i  \  2  .-2X201  ♦  0a>.  ^  .-axa^i  -  0»> 

[a  a 


»ja*ZL3g«-a»*i  +j£*ili.-o0, 

•ji  V  2 


+  *yl  "y2  .-2X0!  -  0a> 

•jl 


)]l 


*W  *yl1.»r«i  ♦/»*> 

4*J> 


where  F(0j)  is  given  by  Eq.  (E7b).  Note  that  >  a  j  since  «j<  x^. 


El.  2  Limiting  Forms  for  Small  Arguments 

In  this  section,  we  evaluate  Z.  (xj,  Xg)  and  2.  (xj,  Xj)  '  (Xg)  arhen 
x1#  Xj  •*  0.  The  expressions  are  applicable  for  0  <  Xj  ^  0. 0001. 

Since 


*n'k 


one  sees  that  when 


then  *rn  •*  0  and.  consequently,  as  given  by  Eq.  (E2a)  and  Eq.  (E2b).  Xj  Xg 
We  distinguish  two  limiting  cases: 

(a)  ¥  »  0,  x  -*  0 

In  this  case,  from  Abramowtts  and  Stegun  (p.  360).  when  x  -  0 
J0(x)  ~  1 
Y0  <x>~5?fnx 
H(J>(x)~  -  j*tnx 
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<E8b> 


(E9a) 

-  0. 

(E  10a) 
(ElOb) 

(ElOc) 


and  therefore 


lim  ^“0^*1*  *2^  * 

*j.yo  xv*i*° 


f  JpOti)  H^W 

l  W 


»  i  -  ltm 
<♦1 


AnOtAp  /l  -<*  ) 
i*(ka  >/l  - 11  ) 


•0,  >»0. 


(BlOd) 


(b)  v  i  0,  Ra  (a)  >  0.  *  <*  0 

In  thla  caaa,  from  AbramowBa  and  St  agon7  (p.  1*0),  Baaaal  functtoae  of  tha 
first  kind*  Neumann  and  Hankal  fUnetiona  can  ba  eapreaaed  In  tarma  ofQamma 
functions  (aaa  Abramowita  and  Stegun. T  p.  158)  aa  follows* 


J,  (x) 


jlil 

r<*+ 1) 


(Ella) 


y„(x)-  -  ^r(a)(f)  * 

H(a2)  (X)- jina)(r)*'. 


(Bilb) 

(Bile) 


Substituting  Eq.  (Ell)  into  Eq.  (BSa),  wa  gat 


Um  Zj,(x.,  x.)  •  lim  — --  /3L)  | 

«,.y.  '  x,.V0  '*  ] 


lim 

>  x-0  r< 

xr*2  ° 


l  /k#o  /m^Yf  /ka  >/l  -  I*  \u] 

~  \  a  /r\kA0^/  J* 


nv 

We  follow  the 


0.  a  M.  (BID 

(D, 


same  procedure  to  evaluate  lim  (x1#  x.)  K*  (x-). 

Xj.Xj^O  * 


Namely* 


T4 


(ft)  ft  «  0,  *  -*  0 

Using  Bq.  (Ell),  «•  hftftft 

Ua  *„<«,.  «,»  ■?'  <V*  -J  •»  U«  (l  -  jj) 

W1  * 


if  lax)  1_ 


4ft 


ks  VI  «  I*  4ftk#oVl  -  ia 
•  M?Ub  In——  -  *  -j7  4b  /  4  / - r 

«j.  V0  **  iJ*  *1  «*»  /,‘l  AU  *“  7l'# 


(E1S) 


k*0  i-1 


(b)  I'M,  *  -*Q 

In  «>i«  mi  mim  Bq.  (Bll)  from  Bq.  (BSs) 


r<y  (grty  • 

Lr(ft  + 1)  . 


lim  Zp(x.,  z.)  (x.)  •  Um 

Wo  W« 

V*/  I  \ k»# /TTp /  J 

*£f£t»  [‘  -  (t^)  ]• 


i 


(B14) 


Utilising  IdooUtlfts 
T(s  ♦  1)  »  ft  ! 

T(ft) « ( v- D! 

Too  _  i 


(ElSft) 


(El  5b) 


(Side) 


r<»  + 1)  * 
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(E19b) 


rm  a7r 

K*(x)~  •"*  ®k  (  '  arg  *!<  ~ 


where 


1  .  HZ  1  +  <«-  -  9)  _  </*  ZJIULZ  9)'(M-28>  +  _ 

8x  2 !  (8x)2  3 !  (8x)3 


(E19c) 


.  !  +  EzA  +  {JLlMJLzlM  ,  -9)(g-25)  + 

^  8x  2!  (8x)2  3!  (8x)3 


(E19d) 


fl  a  4v  . 


(fi!9e) 


Using  Eq.  (E19).  we  have 


Kv<v]  K»<v 


■  —  <s12  *k2  -  '  Xl)  <EJ°> 

•*»  “kJ 


where  again  subscripts  1  and  2  are  associated  with  arguments  Xj  and  Xg#  respec¬ 
tively. 

E2.2  Uniform  Asymptotic  Expansions  for  Large  Orders 
From  Eqs.  (9. 7. 7)  and  (9. 7.  8)7 


Iy  (x)  ~ 


1  e» "I 


s/2nv  (1  +x2)1/4  8i 


(E21a) 


Ky(x)~  jZL 
v  2v 


a"  «I  +  xV/4^ 


(E21b) 


where 


L  Uk<t' 
1  k=l  *'K 


(E21c) 
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^  =  1  +  z  (-l)k  .  <®m> 

k=l 

When  v  -*■  +  «,  these  expansions  hold  uniformly  with  respect  to  x  in  the  sector 
jarg  x|  ^  tt/2  -  £ ,  where  *  is  an  arbitrary  positive  number.  Here 


t  = 


-4 - 


"  ■  Jl  * ?  +*n  i'V^.x* 

and  u^ft)  are  given  by  Eq,  (E5), 

Using  Eq.  (E21)  to  Eq.  (E22),  we  can  write 


i*frl>  1  11 

I„(x2) - -  Kv  (xg)  K„<x?)  =  —  === 

KvtXj)  J  21/^1+xJ 


2 


.  3U  Sk2  e-r2i/072  -  »?i) 
^1 


where 


»2  n  +  Tl  +  x^) 


"2  -  "l  =  ^  +x2  -  ^  +  *1  +4»  , - -  . 

»,(tt  *lt^  I 


5*2. 3  Limiting  Forms  for  Small  Arguments 
Here  we  evaluate 


(S22a) 


<E22b) 


lim 

*1*V°  L 


M*l>  I 


where,  again,  we  distinguish  two  cases: 


(a)  v  »  0,  x  -*■  0 


Using  Eqs.  (9. 6. 7)  and  (9. 6. 8) 


Iq(x)  ~  1 


KQ(x)  ~  -An  x 


one  sees  that 


f  An  kPo  Jt2  -  1  rz - 

ov  1  -  ■  ■—  jmi  i  (“in  kp0  V8Z  “  D 
hl  l  inka  Vi2  -  1  0 


ka  AnkPn  -  1  *a 

-  _j jn -  lim - u_  ..j|n — 

kPQ  An  ka  v$z  -  1  kPr 


(b)  u  t  0.  x  0 


In  this  case,  from  Eqs.  (9.6.7)  and  (9.6.9) 


I*,(x) ~ 


n  v  + 1) 


1  /X 

w-2(j) rM 


and  consequently 


I|/(xi) 

I,,(I2I  '77",  W  Kuix^ 


Xj.Xg-^Q 


v  . v  -  v 


ta  Eiil  •  ixlQl - J  — r<*>  (-S) 

vo[T<^>  ru  +  »)(fr"J  2  '  2 ' 


a  lim 


fi  -(— )2l>(— )" 2"  | 
2  r<“+‘»x1,V0  l  *  »  J 


(E24a) 


(E24b) 


(E24c)  V 


■  1  - 


(E25a)  3 

•  'si 


(E25b) 


i, {.:•$&$£ 

v,  »..v 


, 

■  j  w! 


79 


Appendix  F 
Element  Pattern  Phase  Center 

. . •  -S.'-x'  '*•- 

'  ■"  :'^?3 

The  far  field  due  to  a  singly  excited  element  in  a  match-terminated  array  en¬ 
vironment  is  c&ea tigpUfe s*>k: •  ' ysr  '5* 


,(e) 

J9 


xfT 

(r,  9 ,  g)  * \l  "" 
v  4«rZ. 


Vinc 


•g*  <•> 


g  *  <g*  ♦ ). 


The  phase  reference  point  is  at  r  *  0,  that  is.  at  the  origin  of  the  cylindrical  co¬ 
ordinate  system  (P,  g,  s)  shown  in  Figures  1  and  2.  It  was  found  that  the  flattest 
element  pattern  phase  in  the  broadside  region  is  obtained  when  the  phase  reference 
point  is  at  (p^  ■  a  +  0.8a,  9  ■  9<f .  g  •  OP )  the*  is.  0.2s  below  the  dipole  center. 
Therefore,  it  is  desirable  to  represent  the  far  field  at  (r,  9 , 4 )  with  respect  to  the 
new  phase  reference  point  ( Ppjj,  9  *  90* ,  g  »  <P ).  '•*’** 


7*-i 


Denoting  this  field  by  (R.  9  ,  g )  where 


R  »  r  -  sin  g  cos  g 


(F2a) 


is  the  distance  between  the  new  phase  reference  point  ( Pph»  9  ■  90P .  4  ■  0* )  and 
the  field  observation  point  (r.  9,4 )•  we  may  write 


.  »  u\  JkP  «i»  0  coe  g 

Bj*>  (R.  g.  (r.  g.g)e  ** 


(F2b) 


•  •:|S 
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and  Eg**  (r,  0 ,  ^ )  i»  given  by  Bq.  (FI).  In  Eq.  <F2b),  we  assumed  that,  so  far  as 
amplitude  of  the  field  is  concerned,  l/ra  1/R. 

Substitution  of  Eq.  (FI)  into  Eq,  (F2b)  yields 


E^e)  <R,  9.  *) 


e.4> ) 


(F3a) 


where 


J«>/n 
*6  '  ph 


-ik  p  ,  sin  9  cos  4 
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(F3b> 


and 
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Element  Pattern  of  an  Axial  Dipole  in  a  Cylindrical 
Phased  Array,  Part  2:  Element  Design 

and  Experiments 


1.  INTRODUCTION 

This  part  describes  the  experimental  effort  and  presents  measured  data  that 
strongly  support  the  validity  of  the  theoretical  model  over  a  meaningfully  large 
range  of  element  and  array  parameters.  These  experiments  are  in  support  of  the¬ 
oretical  work  presented  previously*  (see  also  Pert  1). 

The  array  construction  is  described  in  Section  2,  followed  by  the  design  and 
matching  of  the  folded  dipole  radiator  in  SeeMo*  3.  Experimental  amplitude  ele¬ 
ment  patterns  are  presented  in  Section  4  with  theoretical  results  superimposed, 
showing  good  agreement  in  spite  of  the  feet  thatthe  theory  addresses  a  simple  di¬ 
pole  and  the  measurements  were  performed  on  a  folded  dipole.  Measurements 
substantiating  the  theoretical  results  obtained  for  th/ijdment  phase  pattern  are 
discussed  in  Section  5.  Array  beam  forming  effects  are  briefly  considered  in 
Section  6  along  with  polarisation  purity.  The  influence  of  dipole  feed  line  interac¬ 
tions  on  the  radiation  pattern  are  pointed  out. 


(Received  for  publication  31  January  1984) 

1.  Herper,  J.  C. ,  Hessel.  A.„  and  Tomasic,  B.  (1985)  Element  pattern  of  an 

axial  dipole  in  a  cylindrical  phased  array.  Part  1:  Theory.  Part  2:  Element 
design  and  experiments,  IEEE  Trans.  Antennas  Propag. ,  AP-33. 


2.  ARRAY  DESCRIPTION 


An  experimental  test  bed  was  fabricated  to  investigate  the  sidelobe  and  asi- 
muth  phase  scan  capabilities  of  cylindrical  array  antennas.  This  array  was  used 
to  measure  the  element  patterns  over  a  frequency  range  that  exceeded  the  compo¬ 
nent  design  bandwidth.  The  array  is  described  in  this  section*  including  the  feed 
network  and  dipole  geometry. 

The  array  consists  of  80  axial  column  networks  covering  a  150*  arc  of  a  cyl¬ 
inder  as  shown  in  Figure  1.  The  diameter  (2a)  of  the  cylinder  at  the  dipole  ground 
plane  is  29  wavelengths  (all  dimensions*  unless  stated  otherwise,  are  referenced 
to  the  center  frequency  fc).  The  column  networks  extend  approximately  15  wave¬ 
lengths,  with  22  axial  dipoles  spaced  at  0. 67  wavelength.  This  corresponds  to  an 
elevation  incipient  endfire  grating  lobe  angle  of  24*  at  the  high  end  of  the  ±  6  per¬ 
cent  operating  band  for  which  the  array  had  been  designed.  The  column  networks 
are  assembled  into  an  array  with  a  locally  rectangular  lattice  having  circumferen¬ 
tial  spacing  of  0. 5  wavelength  at  the  high  end  of  the  frequency  band.  This  spacing 
was  chosen  to  provide  adequate  suppression  of  the  wide  angle  circumferential  grat- 


Figure  1.  Cylindrical  Array  Test  Bed  Antenna 
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tag  lobes  characteristic  of  fully  excited  cylindrical  arc  arrays  with  large  aperture 
utilisation. 2 3 

The  column  networks  are  microwave  printed  circuits  of  two  types.  The  first 

« 

type  consists  of  22  axially  polarised  folded  dipole  radiating  elements,  each  fed  by 
a  half -wavelength  differential  line  length  balun  and  matched  to  its  own  input  line  by 
a  multiple  step  transformer.  These  radiators  are  used  to  measure  patterns  in  the 
conventional  manner  on  a  far  field  range.  The  other  type  of  column  network  incor¬ 
porates  a  microwave  printed  circuit  that  generates  an  axial  illumination  taper  func¬ 
tion  for  low  sidelobe  patterns.  This  type  of  network,  originally  designed  for  the 
array  studies,  is  useful  for  measuring  circumferential  element  patterns  through 
the  -dal  beam  peak  by  virtue  of  its  increased  isolation  from  ground  reflections  on 
the  measurement  range.  Also,  increased  line  source  gain  over  that  of  the  single 
dipole  improves  the  measured  signal -to-noise  ratio  permitting  circumferential 
patterns  to  be  measured  to  large  angles  before  becoming  noise  limited.  The  meas¬ 
ured  noise  level  of  the  single  dipole  element  patterns  is  approximately  lOdB  poorer 
than  the  noise  of  the  column  patterns,  in  good  agreement  with  the  gain  improve¬ 
ment.  Except  for  the  noise  level,  measurements  with  both  types  of  networks  for 
the  same  conditions  provided  identical  results. 

The  column  networks  utilise  the  folded  dipoles  and  baluns  discussed  above, 
fed  by  a  22:1  corporate  power  divider  designed  to  produce  a  -35 dB  Taylor  illumi¬ 
nation  using  compensated  Wilkinson  couplers.  The  symmetric  and  four -port  re¬ 
sistive  nature  of  these  coupler  makes  them  particularly  useful.  They  provide  a 
uniform  power  split  in  amplitude  and  phase  as  a  function  of  frequency,  excellent 
wideband  match,  and  high  isolation.  Each  of  the  dipoles  and  the  column  network 
input  is  effectively  isolated  frcm  the  mismatches  created  at  other  dipoles  during 
variations  of  scan  and  frequency.  Asa  result,  any  dipole  in  the  array  sees  an  ex¬ 
cellent  feed  match  over  all  scan  and  frequency  conditions.  The  printed  circuits 
are  copper-clad  Kapton,  sandwiched  between  aluminum  ground  planes  using  low  di¬ 
electric  constant  foam  spacers.  The  response  of  the  couplers  and  transmission 
lin*s  is  flat  in  amplitude  and  phase  over  the  12  percent  operating  band.  Measure¬ 
ments  on  the  80  column  networks  over  the  operating  band  showed  rms  phase  and 
amplitude  tolerances  of  2.1s  and  0. 15dB,  respectively.  A  picture  of  the  network 
is  shown  in  Figure  2  with  the  protective  foam  partially  removed  to  display  the  di¬ 
poles. 


2.  Hessel,  A.  (1972)  Mutual  coupling  effects  in  circular  arrays  on  cylindrical 

surfaces  -  aperture  design  implications  and  analysis.  Phased  Array  Anten¬ 
nas.  Artech  House,  Dedham,  Mass. ,  p.  273. 

3,  Jasik,  H.  (1959)  Antenna  Engineering  Handbook,  McGraw-Hill,  New  York, 

p.  3-13  to  3-15: -  - 
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Figure  2.  Column  Network  With 
Protective  Foam  Partially  Re¬ 
moved  Showing  Dipole  Radiator  a 


3.  DIPOLE  DESIGN 

A  folded  dipole  fed  by  a  differential  line  length  balun  was  chosen  because  of  its 
extensive  use  in  planar  arrays  and  good  correlation  to  theoretical  performance  for 
those  applications.  A  previously  used  design  was  scaled  to  the  desired  operating 
band,  and  the  element  was  matched  in  a  planar  waveguide  simulator  using  both  re¬ 
flection  and  transmission  measurements.  An  outline  diagram  of  the  radiator  is 
shown  in  Figure  3.  The  dipole  is  0. 42  wavelength  long  and  is  located  at  0.24 
wavelength  above  the  cylindrical  ground  at  center  frequency.  It  is  fed  by  a  two- 
wire  line  in  the  free  space  region  and  matched  to  two  stripline  conductors  of  the 
same  impedance.  The  balun  is  formed  in  this  stripline  by  adding  a  half-wave- 
length  of  line  on  one  feed.  The  two  lines  are  then  reactively  combined  into  a  70J3 
characteristic  impedance  stripline  used  for  the  remainder  of  the  printed  circuitry. 

3 

As  described  by  Jasik,  the  folded  dipole  performs  in  a  fashion  similar  to  the 
ordinary  center-fed  dipole  with  the  folded  lines  providing  an  impedance  trans¬ 
former  between  the  dipole  and  its  feed  line.  The  nearest  to  broadside  single  - 
element,  single-mode  45*  H -plane  simulator  shown  in  Figure  4a  was  employed 
for  matching  the  dipole.  The  simulator  was  excited  at  the  Mid  opposite  the  dipole 
with  an  orthogonal  mode  transducer.  This  provided  a  conveniert  means  to  meas¬ 
ure  and/or  terminate  the  cross-polarised  radiation.  Also  shown  in  Figure  4b  is  a 
test  dipole  in  its  mating  ground  plane  fixture.  Figure  4c  is  a  diagram  of  the  ele¬ 
ment  lattice.  For  this  configuration,  the  simulator  height  is  the  same  as  the 
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Figure  3.  Outline  Diagram  of 
Folded  Dipole  Radiator  With 
Stripline  Feed 


axial  dipole  spacing  and  the  simulator  width  is  50  percent  greater  than  the  circum¬ 
ferential  dipole  spacing.  The  dipole  is  located  one**hlrd  the  simulator  width  from 
the  left  side  wall.  The  simulator  is  exact  for  an  element  in  a  planar  array.  It  ap¬ 
proximates  the  active  dipole  impedance  at  the  average  scan  in  the  arc  excited  for 
broadside  beam  formation  in  a  large  cylindrical  array.  The  resulting  impedance 
match  over  the*  6  percent  design  operating  band  is  shown  in  Figure  5.  A  voltage 
standing-wave  ratio  (VSWR)  of  better  than  1. 3:1  was  achieved  over  the  operating 
band. 

4k  MEASURED  AMPLITUDE  PATTERNS 

Dipole  element  amplitude  patterns  were  measured  cm  a  87 5 -ft  antenna  range. 
Extensive  measurements  were  performed  in  the  ±  6  percent  design  operating  band. 
Selected  measurements  were  also  made  up  to  30  percent  above  the  design  operat¬ 
ing  frequency.  All  measured  dipole  patterns  correlated  very  well  with  the  theo¬ 
retically  predicted  performance. 

The  experimental  data  required  somewhat  modified  computations  compared  to 
those  presented.  First,  since  the  experimental  dipole  is  matched  ai  OP  in  the 
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Figure  4.  Waveguide  Simulator  for  45*  H -Plane,  (f  SE -Plane  Scan, 
(a)  Simulator  body  with  orthogonal  moade  tranaduaer.  (b)  Teat 
dipole  in  fixture,  (c)  Relationship  of  array  lattice  to  simulator 
size  and  location 


E-plane  and  45*  in  the  H -plane,  the  theoretical  data  shown  here  were  recomputed 
to  reflect  this  condition.  In  general,  this  effect  causes  a  reflection  mismatch  of 
about  0. 5  dB  at  broadside,  which  flattens  the  central  portion  of  the  pattern  out  to 
the  match  angle  and  slightly  narrows  the  patterns  for  angles  greater  than  the 
match  angle.  The  second  change  incorporated  into  the  computer  simulation  arises 
from  the  fact  that,  in  the  experiments,  the  physical  dimensions  are  held  constant 
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Figure  5.  Measured  Dipole  Impedance  in 
Waveguide  Simulator 


-while  the  frequency  varies.  This  leads  to  a  set  of  normalised  element  and  array 
dimensions  that  increase  linearly  with  the  impressed  source  frequency.  Tables 
are  provided  on  each  pattern  indicating  the  normalised  dimensions  of  the  antenna 
at  the  measurement  frequencies.  The  symbols  are  defined  in  a  list  at  the  end  of 
this  report. 

The  measured  H -plane  pattern  at  the  center  frequency  (f  *  fc)  is  shown  as  the 
solid  curve  in  Figure  6.  This  condition  corresponds  to  a  circumferential  spacing 
between  dipoles  of  0. 48  wavelength.  The  curve  shows  the  element  gain  in  dB  nor¬ 
malised  to  that  of  the  unit  cell  area  versus  the  circumferential  field  angle  d  mess 
ured  from  the  element  broadside  as  the  abscissa.  This  is  a  principal  H -plane  col 
umn  pattern  (taken  normal  to  the  cylinder  axis*  that  is.  at  $  •  9<f  with  respect  to 
the  cylinder  axis).  Measurements  were  terminated  at  -40  dB  below  the  beam  peak 
( #  «  ±  12QT )  because  of  test  and  sits  instrument  at  ion  noise  limitations.  Superim¬ 
posed  on  this  plot  are  theoretical  points  represented  by  dots.  The  excellent  agree 
meut  is  evident;  this,  in  spits  of  the  tact  that  the  theoretical  model  employs  a 


simple  strip  dipole  of  uealigibie  redial  thickness,  while  the  actual  folded  dipole  is 


printed  on  radial  boards.  The  length  L  of  the  theoretical  model  was  assumed 
equal  to  that  of  the  top  arm  of  the  folded  dipole. 

A  similar  set  of  data  at  a  frequency  4  percent  higher  showing  excellent  corre¬ 
lation  between  experiment  and  theory  is  shown  in  Figure  7.  This  result  is  of  in- 


CIRCUMFERENTI AL  ANGLE  <t>.  (DEG) 


Figure  6.  Measured  and  Theoretical  Amplitude 
Patterns  at  Center  Frequency 


CIRCUMFERENTIAL  ANGLE  0.  (DEG) 


Figure  7.  Amplitude  Patterns  4  Percent  Above 
Center  Frequency 
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terest  for  comparison  to  slot  data  provided  in  the  literature.  This  frequency 
provides  an  effective  element  spacing  of  0.  50  wavelength.  In  contrast  to  a  slot  ar¬ 
ray  with  this  spacing,  the  dipole  patterns  show  no  discernible  element  pattern 
ripples  in  the  forward  region.  For  further  comments  in  this  connection,  see 
Part  1  of  this  report. 

Data  for  a  circumferential  element  spacing  of  0.  594  wavelength  at  f  =  1, 2fc, 
is  shown  in  Figure  8.  To  accentuate  the  pattern  ripple  in  the  lit  region,  this  graph 
is  plotted  with  the  ordinate  on  a  linear  (voltage)  scale,  again  normalized  to  the 
unit  call  gain.  Correlation  of  the  experimental  data  (solid  curve)  witn  the  theory 
(dots)  is  again  remarkably  good,  especially  since  the  test  frequency  is  more  than 
20  percent  above  the  center  frequency  at  which  the  components  were  designed  to 
operate. 

A  series  of  measurements  presented  in  Figure  9  shows  the  dipole  pattern  per¬ 
formance  parametrically  at  a  number  of  conical  cut  angles  relative  to  the  cylinder 
axis.  These  patterns  are  for  polar  angles  Q  of  90“  (Figure  9a),  80°  (Figure  9b), 
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CIRCUMFERENTIAL  ANGLE  <t>,  (DEG) 


Figure  8.  Amplitude  Patterns  22.6  Percent 
Above  Center  Frequency 


4.  Sureau,  J.  C.,  and  Hessel,  A.  (1971)  Element  pattern  for  circular  arrays  of 
waveguide-fed  axial  slits  on  large  conducting  cylinders,  IEEE  Trans.  An¬ 
tennas  Propag.  AP-19:64-76. 


91 


CIRCUMFERENTIAL  ANGLE  <l>,  (OEG) 
(a) 


CIRCUMFERENTIAL  ANGLE  0,  (DEG) 
(b) 


Figure  9.  Amplitude  Patterns  at  Various 
Elevation  Angles,  (a)  9  =  0*  •  (b>  $  *  1  OP 


70"  (Figure  9c),  and  60"  (Figure  9d).  In  contrast  to  the  previous  patterns,  which 
were  column  patterns,  these  element  patterns  are  of  lower  absolute  gain  by  about 
lOt'B  and  are  thus  affected  by  measurement  system  noise  at  a  10 dB  higher  level. 
Pejk  gain  decreases  with  elevation  angle  at  nearly  a  cosine  voltage  pattern  as  ex¬ 
pected  for  an  E-plsne  dipole  scan  when  both  the  unit  cell  directivity  and  scan  mls- 
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CIRCUMFERENTIAL  ANGLE  <#>,  (DEG) 
(c) 


CIRCUMFERENTIAL  ANGLE  O.  IOEGI 
(d) 


Figure  9  (cont. ).  Amplitude  Patterns  at  Various 
Elevation  Angles,  (c)  0  *  2(f  .  (d)  §  *  30P 


match  gain  loss  are  included.  The  patterns  show  a.  general  circumferential  broad¬ 
ening  with  increasing  elevation  scan.  As  indicated  by  the  superimposed  dots,  the 
actual  performance  is  very  close  to  the  theoretically  predicted  performance. 
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5.  PHASE  PERFORMANCE 


In  Part  1  of  this  report,  the  phase  response  of  the  dipole  was  found  to  be  flat 
to  within  a  few  degrees  ovei  a  large  portion  of  the  lit  region  when  the  phase  refer¬ 
ence  distance  from  the  ground  plane  is  taken  about  80  percent  of  the  distance  from 
the  ground  to  the  dipole  strip.  Instead  of  measuring  the  element  phase  directly,  a 
rather  difficult  undertaking  on  a  far-field  range,  an  indirect  focusing  measure¬ 
ment  was  performed. 

The  focusing  experiment  consists  of  the  following  tasks.  A  convenient  illumi¬ 
nation  function  is  impressed  on  the  array.  Far-field  patterns  are  measured  as 
the  element  phase  excitation  is  varied  by  incrementing  the  apparent  radius  of  the 
cylindrical  array  illumination.  This  simulates  the  variation  of  the  dipole  phase 
center  location.  An  error  in  the  phase  center  location  impresses  a  quadratic 
phase  error  on  the  array  collimation  which,  in  turn,  manifests  itself  as  pattern 
defocusing.  From  a  practical  viewpoint,  approximately  the  same  results  may  be 
obtained  more  easily  by  varying  the  test  frequency  without  readjusting  the  array 
collimation.  This  effectively  simulates  the  desired  collimation  variation. 

The  experiment  begins  with  the  radius  set  to  the  actual  dipole  location.  As 
the  radius  (equivalent  dipole  radial  distance)  is  reduced,  the  nulls  of  the  close-in 
sidelobes  increase  in  depth,  indicating  an  improvement  in  the  array  focusing. 

For  a  radius  corresponding  to  70  percent  of  the  distance  from  the  ground  plane  to 
the  dipole,  these  nulls  reach  maximum  depth.  Null  depths  at  this  point  are  about 
lodB  better  than  at  the  beginning  of  the  experiment.  Further  reduction  of  the  ra-- 
dius  only  degrades  the  null  depths,  confirming  that  the  optimum  location  has  been 
exceeded.  These  measured  results  are  shown  graphically  in  Figure  10. 

The  theoretically  predicted  phase  center  location  thus  agrees  very  well  with 
the  results  of  the  focusing  measurement.  This  good  focus  held  over  a±  6  percent 
test  band  and  with  circumferential  phase  s  an  of  the  array.  The  agreement  with 
theory  is  further  accentuated  since  the  folded  dipole  feed  employed  in  the  meas¬ 
urements  can  be  expected  to  have  a  slightly  different  phase  center  location  than 
the  flat  strip  dipole  investigated  in  the  theory. 


6.  ARRAY  BEAM  FORMING  CONSIDERATIONS 

An  important  effect  in  cylindrical  arrays  is  the  formation  of  residual  delete¬ 
rious  grating  lobes.  This  effect  is  caused  by  the  edge  subarrays  of  the  aperture, 
which  need  to  be  given  large  progressive  phase  delays  to  collimate  the  beam  and 
which  have  large  enough  spacing  to  generate  a  grating  lobe.  Using  the  precise 
element  patterns  generated  by  the  methods  described  in  Part  1,  the  grating  lobes 
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Figure  10.  Location  of  Dipole  Phase  Center 


can  be  accurately  predicted.  These  lobas  generally  occur  at  azimuth  angles 
greater  than  90*  from  broadside,  may  have  large  amplitude  compared  with  the  de¬ 
sired  sidelobe  levels,  and  are  usually  very  broad  in  angular  extent.  These  effects 
are  demonstrated  in  the  literature.  * '  *  The  oii^  design  tool  effective  in  con¬ 
trolling  these  grating  lob  as  is  the  reduction  of  the  circumferential  spacing  between 
the  elements  of  the  array.  Figure  11  shows  the  result  of  a  pav  ametric  study  of 
this  effect.  A  cylinder  of  constant  radius  (ka  -  140)  is  considered,  and  the  grating 
lobe  level  is  plotted  versus  circumferential  element  spacing.  A  120*  sector  of  the 
array  is  assumed  to  be  excited  with  a  low  sidelobe  illumination  at  broadside  and  * 
15*  circumferential  phase  scan,  which  corresponds  to  75°  scan  for  the  edge  ele¬ 
ment.  When  the  beam  is  broadside  to  the  excited  aperture,  the  far  out  sidelobes 
computed  in  the  region  where  the  grating  lobe  appears  approached  -70dB  for  close 
element  spacings.  For  element  spacing  near  0.  50  wavelength,  the  residual  grat¬ 
ing  lobe  level  generated  is  negligible.  As  the  element  spacing  approaches  0. 60 
wavelength,  lobes  well  above  -40  dB  (relative  to  the  main  beam)  are  produced. 

For  many  applications,  such  grating  lobe  levels  would  be  unacceptable. 

Another  deleterious  mutual  coupling  and  curvature  effect,  discussed  in  Part  1 
of  this  report,  results  from  the  behavior  of  the  element  pattern  as  a  function  of 
element  spacing.  At  element  circumferential  spacings  corresponding  to  0.  5  wave- 


5.  Sureau,  J.C.,  and  Hessel,  A.  (1972)  Realized  gain  function  for  a  cylindrical 

array  of  open-ended  waveguides.  Phased  Array  Antennas.  Artech  House, 
Dedham,  Mass.,  p.  315. 

6.  Provencher,  J.H.  (1972)  Conformal  arrays  on  surfaces  with  rotational  sym¬ 

metry,  Phased  Array  Antennas,  Artech  House,  Dedham,  Mass.,  p.  301. 
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Figure  11.  Grating  Lobe  Level  of  a  Low  Sidelobe 
Pattern  Excited  by  a  120*  Cylindrical  Array  Arc 


length,  a  smooth  dipole  pattern  is  observed.  For  spacings  greater  than  this,  a 
ripple,  due  to  creeping  waves  propagating  around  the  cylinder,  appears  in  the 
broadside  region  of  the  element  pattern. 2  In  a  cylindrical  array,  this  ripple  can 
act  like  a  periodic  error  across  the  array,  and,  at  some  frequencies,  may  lead  to 
the  generation  of  anomalies  in  the  near -in  array  sidelobes.  Indeed,  element  rip- 
pit  levels  of  a  few  tenths  of  a  dB  peak-to-peak.  could  cause  sidelobe  peaks  on  the 
order  of  -40  dB  below  the  array  beam  peak.  *  nis  effect  on  the  sidelobe  level  is 
generally  difficult  to  predict  because  it  depends  on  the  proper  correlation  of  the 
periodicity  of  the  element  pattern  ripple  with  the  effective  element  spacing,  each 
of  which  depends  on  frequency.  Many  frequency  patterns  must  be  evaluated  to 
guarantee  suppression  of  this  effect  to  acceptable  levels.  If  measured  patterns 
.are  employed  to  predict  array  performance,  care  must  be  exercised  in  separating 
the  actual  ripples  obtained  on  a  full  cylindrical  array. from  those  caused  by  exper¬ 
imental  anomalies  such  as  range  multipath  or  edge  effects  due  to  finite  subarray 
test  units. 

It  is  also  useful  to  compare  the  dipole  performance  with  that  of  a  slot  radia¬ 
tor.  4  For  equivalent  size  cylinders  and  element  spacings,  the  dipole  has  consid¬ 
erably  less  pattern  ripple  than  the  slot.  This  is  attributed  to  the  imaging  property 
of  the  dipole  above  its  ground  (see  Part  1  of  this  report).  Also,  the  dipole  has  a 
slightly  broader  pattern  due  to  the  weaker  grating  lobe  effect  as  compared  to  the 
slot.  These  features  may  favor  the  dipole  in  array  designs  that  utilize  significant 
portions  of  a  cylindrical  arc  or  where  low  sidelobes  are  desired. 

As  discussed  in  Part  1  of  this  report,  the  longitudinal  current  approximation 
for  the  axial  strip  dipole  does  not  predict  any  circumferential  far  E -field.  Array 
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pattern  measurements  confirm  this  analysis.  The  crosB-polarization  was  meas¬ 
ured  over  a  large  region  about  the  beam  in  the  circumferential,  axial,  and  skew 
planes.  The  cross -polarization  level  was  everywhere  better  than  10  dB  below  the 
corresponding  co -polarization  level  at  the  same  angle  and  was  generally  as  good 
as  the  cross -polarized  level  of  the  site  transmit  dish.  Thus,  no  experimental  evi¬ 
dence  exists  for  coupling  to  the  cross-polarized  radiation  within  the  measurement 
region.  This  indicates  that  other  structures  such  as  the  feed  lines  and  the  folded 
feed  that  could  couple  to  the  cross-pclarized  radiation  generate  no  significant  cir¬ 
cumferential  electric  field. 

It  has  been  noted  in  the  literature  that,  under  certain  circumstances,  the  di¬ 
pole  feed  lines  can  resonate  with  the  dipole  to  form  pattern  blind  spots  in  the  E- 
7  8 

plane.  '  Data  and  discussion  of  this  effect  are  scanty,  difficult  to  obtain,  and 
often  limited  to  verbal  presentations.  Pattern  data  in  the  axial  (E)  plane  were 
limited  to  30°  because  of  turntable  loading  considerations.  No  resonance  effects 
were  encountered  out  to  this  point. 


7.  CONCLUSION 

The  measured  results  for  the  element  pattern  of  an  axial  dipole  element  in  a 
cylindrical  phased  array  lead  us  to  the  following  conclusions: 

1 .  The  modal  analysis  developed  in  Part  1  of  this  report  agrees  very  well 
with  the  measured  data  indicating  that  the  method  and  the  numerical  results  are 
very  accurate. 

2.  The  assumption  that  the  feed  lines  have  little  effect  is  justified.  For  the 
scans  and  frequencies  investigated,  no  differences  were  discernible  between  the 
theory  (no  feed  modeled)  and  the  experiments  with  the  folded  end -fed  dipole  and 
two -wire -line  feed. 

3.  A  high  polarization  purity  and  wide  bandwidth  dipole  and  feed  are  achiev¬ 
able. 

4.  The  accurate  knowledge  of  element  pattern  amplitude  and  phase  afforded 
by  this  method  makes  possible  the  precise  prediction  of  cylindrical  array  perfor¬ 
mance  over  a  full  360P  azimuth  range. 
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1 


Nomenclature 


e 

* 

$ 

#o»  ^  o 

ka 

b 

d 

L 

s 

vr 

h 

N 

f 

fc 


Polar  angle  measured  from  the  cylinder  axis 
Circumferential  angle  measured  from  element  broadside 
=  (  90*  -  0  )  -  Elevation  angle  measured  from  normal  to  axis 
Beam  pointing  direction  for  active  phased  array 
Array  circumference 
Circumferential  element  spacing 
Axial  element  spacing 
Dipole  length 

Dipole  spacing  from  cylindrical  ground 
Dipole  width 
Dipole  feed  gap 

Number  of  dipoles  in  each  ring  (192) 

Operating  frequency 
Design  center  frequency 


NOTE:  All  lengths  in  wavelengths  at  the  test  frequency. 
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